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TECHNICAL NOTE 2050 


A METHOD FOR CALCULATING THE SUBSONIC STEADY-STATE 
LOADING ON AN AIRPLANE WITH A WING OF ARBITRARY 
PLAN FORM AND STIFFNESS. 


By W. L. Gray and K. M. Schenk 
SUMMARY 


A method for computing the steady-state span load distribution on 
an elastic airplane wing for specified airplane weights and load factors 
is given. The method is based on a modification of the Weissinger 
L-method and applies at subcritical Mach numbers. It includes the 
effects of external stores and fuselage on the spanwise loading. Modi- 
fications are outlined for treating tail-boom and tailless airplane con- 
figurations and for calculating the divergence dynamic pressure of a 
swept wing with a large external store. A method is also outlined for 
reducing wind-tunnel data to obtain effective aerodynamic coefficients 
which are free of model flexibility effects. The effects of Mach num- 
ber can readily be evaluated from the aerodynamic coefficients thus 
obtained. 


INTRODUCTION 


The inclusion of the effects of flexibility in the solution of the 
spanwise airload distribution applied to a wing of arbitrary plan form 
and stiffness distribution has increased the complexity of analysis over 
that for a rigid wing. The methods that are available at the present 
time are generally concerned with the calculation of loading on an iso- 
lated flexible wing rather than the more practical case not only where 
the effects of fuselage and nacelles on the spanwise loading must be 
taken into account but also where the total lift on each of the major 
components must be considered simultaneously in order to determine the 
wing loading at a specified load factor. A method for including such 
effects without recourse to iterative procedures for steady-state flight 
conditions and subcritical Mach numbers is presented in this paper. The 


lBased on Boeing Airplane Company Document No. D-1062h, "A Matrix 
Solution for the Subsonic Steady State Aeroelastic Loading on Airplanes" 
by W. L. Gray and K. M. Schenk, June 1, 1951. Acknowledgement is made 
to Messrs. Paul W. Harper and John B. Garvin of the NACA for extensive 
work in editing and revising this document. 
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equations are derived so that the spanwise airload distribution can be 
expressed in matrix form in terms of influence coefficients for aero- 
dynamic induction and structural deflection in a manner similar to that 
employed in reference l. 


The basic method is outlined in the body of the paper. Included 
in appendixes are details of the various derivations, the expansion of 
the basic equations to include fuselage interference and store load 
effects, the modifications for tail-boom and tailless configurations, 

a method for determining divergence dynamic pressures for swept wings 
with large external stores, a method for reducing wind-tunnel data to 
obtain effective aerodynamic coefficients which are free of model Plexi- 
bility effects, and a method for obtaining compressibility corrections. 


SYMBOLS 


The following symbols appear in the body of this report. Addi- 
tional symbols which appear only in the appendixes are defined as they 
are introduced. 


b wing span, in. 

Cp o fuselage lift coefficient in presence of wing at ap = 0, 
Lp/as 

mo fuselage pitching-moment coefficient about C/4 point in 
presence of wing at a, = 0, Mp/qs¢e 

(erg). rate of change of fuselage lift coefficient with a, 
per radian 

(Cin) rate of change of fuselage pitching-moment coefficient 

Co with ap, per radian 
c local chord parallel to plane of symmetry, in. 
c wing mean geometric chord, in. 


Cm. wing section pitching-moment coefficient 
EI 


effective value of product of modulus of elasticity and 
wing section beam bending moment of inertia, lb-in.2 


T 


effective value of product of shear modulus of elasticity 
and wing section polar moment of inertia, lb-in.@ 
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semispan of horseshoe vortex, in. 
airplane rolling moment of inertia, lb-in.? 


airplane pitching moment of inertia, lb-in.? 


fuselage lift in presence of wing, [ereo + (Ct) ar |as, lb 
a 


Wing running lift per inch of span perpendicular to plane 
of symmetry, lb/in. 


fuselage pitching moment about 2/4 point in presence of 
+ Sc, in-1b 
wing, Gne, * (eme) er 956, 


airplane rolling moment, in-1b 


elastic-airplane rolling moment caused by unit aileron 
deflection, in-1b 


elastic-airplane damping moment in roll caused by unit 
wing-tip helix angle (pb/2V),, in-1b 


elastic-airplane rolling moment caused by unit rolling 


angular acceleration acting on the wing distributed 
inertia, in-1b 


elastic-airplane rolling moment caused by unit spoiler 
deflection, in-1b 


two-dimensional lift-curve slope per radian, including com- 
pressibility effects, for sections parallel to plane of 
symmetry 


airplane load factor, positive when inertia loads are downward 


balancing tail load, positive upward, lb 


airplane rolling velocity, radians/sec 

airplane rolling angular acceleration, radians/sec® 
dynamic pressure, lb/sq in., pv*/2 

radial distance from vortex core, in. 


wing area, sq in. 
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true free-stream velocity, in./sec 
airplane gross weight, lb 


wash velocity induced by line vortex at perpendicular dis- 
tance r from vortex line, positive for downwash, in./sec 


downwash angle at three-quarter-chord point induced by 
vortex system representing wing and its spanwise lift 
distribution 


streamwise distance from pitch reference axis to bound 
portion of horseshoe vortex, positive when vortex is 
to rear of pitch reference axis, in. 


streamwise distance from pitch reference axis to T/4 line, 
positive when C/4 line is to rear of pitch reference 
axis, in. 


streamwise distance from pitch reference axis to airplane 
center of gravity, positive when center of gravity is to 
rear of pitch reference axis, in. 


streamwise distance from pitch reference axis to center of 
pressure of balancing tail load, positive when center of 
pressure is to rear of pitch reference axis, in. 


lateral distance from wing center line, in. 

final angle of attack of section zero-lift line with respect 
to local free-stream direction, a, + Gg + Gg, radians 
(see fig. 1) 

change in section angle of attack due to aerodynamic twists 


and due to all structural twists associated with a flexi- 
ble wing which are not accounted for by the a, term, 


radians (see fig. 1) 


angle of attack of root-section zero-lift line, radians 
(see fig. 1) 


change in section angie of attack due to wing lift distri- 
bution acting on a flexible wing («s = 0 fora rigid wing), 


radians (see fig. 1) 
strength of line vortex, in.2/sec 


aileron deflection, radians 


Fota 
P wer os uu à e ` - 
- ` 5a -. r 
- ^* —e = Y^ * -€-*1 
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Og spoiler deflection 

y dimensionless spanwise station, JE 

9 airplane pitching angular acceleration, positive for nose up, 
radians/sec* 

A. local sweep angle of elastic axis, radians 

Aw equivalent local sweep angie including compressibility 


effects, radians 


p mass density of ambient atmosphere, slugs/cu in. 


(p = 0.114679 x 10-6 1b-sec®/in.4 at standard sea- 
level conditions) 


Matrix notation: 


[ | square matrix, elements of which are designated by.use of 
subseripts; for example, element 813 is in ith row and 
jth colum 


D E row matrix 
i : column matrix 
O 


diagonal matrix, which is a square matrix in which all 
| | elements are zeros except those on the principal 


diagonal 814, oo, zz) * * * Onn 


El aerodynamic-induction or downwash matrix in which ele- 
ments 84 j relate downwash angle at station i to 


unit running lift at station j on wing 


[s2] elasticity matrix in which elements ajj relate changes 
in streamwise angle of attack at station i to unit 
running lift at station j on wing 


[Ss | fuselage image-vortex matrix relating image downwash 
effects at station control points to unit running lifts 
(see appendix E) 


O 
[So | fuselage "overvelocity" matrix (see appendix E) 


-—- - -ai o4 AR o0 ~ vc 7, 
‘ e 
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[1] identity matrix; that is, diagonal matrix in which 
diagonal elements are equal to unity 


PRESENTATION OF METHOD 


In this section of the report the basic equations necessary to the 
method are outlined and discussed in a general wey. Details of the 
derivations are contained in the various appendixes. 


Assumptions 


In the development of the method certain assumptions that are 
common to airfoil theory apply, namely: 


(1) The flow is potential; that is, boundary-layer effects, separa- 
tion, and compressibility shocks are absent or negligible. 


(2) The wing thickness is small. 
(3) A stagnation point exists at the wing trailing edge. 


(4) The angles of attack a are small so that tana * sin a a 
(where a is measured in radians) and cosa 9 1. 


(5) All drag-load effects except those due to nacelles and stores 
are neglected entirely in determining the deformations of the wing used 
in obtaining the equilibrium spanwise airload distribution. 


With regard to the structure the following assumptions are made: 


(1) Camber changes arising from twisting and bending of the wing 
are neglected entirely. 


(2) The elastic twist of the control surface is the same as that of 
the adjoining wing structure. 


(5) The angles of structural deflection 6 are small so that 
tan 6 = sin 6 * 0 (where 8 is. measured in radians) and cos 0 = 1. 


(4) Although the &ngle-of-attack changes, including those due to 
bending and torsional deformations of the wing, are accounted for in 
the determination of the equilibrium spanwise airload distribution on 
the wing, this final airload distribution is applied to the geometry of 
the undeflected wing in computing the bending and torsional moments. 
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Basic Equations 


symmetrical flight conditions.- The fundamental problem involved is 
the development of a series of equations which relate the spanwise lift 


distribution for an arbitrary wing plan form in a given flight condition 
to the properties and attitudes of the individual sections that form the 
wing. 


If the two-dimensional wing is considered first, the following 
relationships for lift and downwash behind an airfoil are available 
from most standard textbooks on aerodynamics: 


l = ovr (1) 
l = Pt p x C (2) 
W. = m (5) 


The circulation [ is taken to be such that, at a specified dis- 
tance r behind the lifting line, the resultant of the downwash veloc- 
ity w, and the flight velocity V is parallel to the section zero-lift 


line; that is, no flow exists normal to the zero-lift line at this point. 
Then, 


y Te (4) 


DP = Mae pe (5) 
Substituting: equation (5) into equation (3) results in 


2 
Wy = To SÅ agy , (6) 
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or 
W: 2 
uo ae (7) 
V ex r 
mo c/2 
In order to satisfy equation (4), the expression Sp in equa- 
T 


tion (T) must be equal to 1.0. Since the theoretical section two- 
dimensional lift-curve slope is equal to 2x, r must equal c/2, which 
is the distance between the lifting line and the three-quarter-chord 
point. 


In the development of the method presented in this report, equa- 
tion (T) is always used in the form 


Um j Ba Ot is 


This simplification requires that the section lift-curve slope mg be 
the two-dimensional value (i.e., the value of the lift-curve slope for 
an unswept two-dimensional wing) and that the location of the downwash 
control point D (see fig. 2) be one-half of the local streamwise chord 
to the rear of the quarter-chord point, or at 3c/%. 


The essential difference between a two-dimensional wing and a wing 
of finite aspect ratio arises from the nonuniform spanwise loading which 
produces the trailing vortices of the finite-aspect-ratio wing. The 
equations presented thus far are considered to apply to the finite- 
aspect-ratio wing when the effects of all the vortices, both bound and 
trailing, have been taken into account. 


Equation (8) in matrix form is 


O 


UM i E " 


This matrix relation represents a series of equations, each applicable 
to a particular station on the semispan of the wing. ‘The values of 


= ie every one of which is affected by the bound and trailing 
5c 
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vortices at all of the wing stations, can be evaluated from 

Ween TREJO 

3c h^ LV 

which, in combination with equation (1), results in 

d = [81] fa} (10) 

V = AR. Yl 

t 3c pr wn 1 


The EN matrix in these equations is the aerodynamic-induction or 
downwash matrix which is derived in &ppendix A. 


Combining equations (9) and (10) gives 


fa... EO - Bley an 
Jue -e en 


The series of equations represented by the matrix equation (12) expresses, 
for any given dynamic pressure, the relationship between the spanwise 


variation of running lift {i}, the final section angles of attack ted , 
and the spanwise variation of the two-dimensional section lift-curve 

O 
slope [Mo |- The effects of wing plan-form geometry are accounted for 


through the elements of the [$1] matrix. The section 1ift-curve slope 


is expressed in the general form mj rather than ex to permit substi- 


tution of actual values when available from scaled-model tests or to 
permit correction for compressibility effects as described in appendix A. 
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The final angle-of-attack variation across the span {ap} can be 


considered to be composed of three essential parts (see fig. 1) 


eepe o 


For a wing free of external stores, the angle of attack fas} caused 


by structural deflection of a flexible wing due to the section lifts 
acting at the section aerodynamic centers is linearly related to the 


matrix Uu by an expression derived in appendix B as 


| fox} = [5] {i} (1) 


The wing geometry and stiffness are accounted for in the structural- 
deflection matrix [S2]- This matrix is based on loadings associated 


with stations which are parallel to the airplane plane of symmetry. In 
a swept wing, however, the structure is usually arranged such that the 
wing boxes are formed between stations approximately perpendicular to 
the elastic axis. In order to obtain a closer representation of the 


loadings and deflections on this actual wing, an EM matrix was also 


derived (see appendix B) and may be substituted for the [S2] matrix in 


equation (14) when desirable. 


The les] matrix of equation (13), as described in detail in 


appendix C, is composed of built-in twist, apparent or aerodynamic twists 
such as those due to interference, control deflection, and angular veloc- 
ities, and all structural twists of an elastic wing which are not accounted 


for in the {ash matrix. 


Although equation (12) is perfectly general, it is not useful in 
the form given for determining the lift distribution on a flexible wing 


since a component of the {ae} matrix is itself a function of the lift. 


If {as is therefore expressed as in equation (1+), equations (12) 


Le E —— — ————— a me A O a aper pe 
A de ^ = 3 m : - » P ~ ce ren ne a a 
me . A » r MM r - z ra ^ - t : 
" E t " " » - 1 * weet -z - P Tode a ^or + F %s k * - + tog? a 
` y i 


y * 
. A . a a pS - A * a > . ` 
- * - * 2.7 a* * e - - ‘ ld a 1 LI - 
` - LES - LED - . . x a . 
na E ML MESS EOS ^, . 45 . nna ~op a ye, " RIS 2 tT n A = 
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and (13) may be combined so as to express the load distribution on a 
fiexible wing in terms of wing root angle of attack and any combination 


of the {og} twists as 


bs] - 63 +6} raO os 


or 


[es E - [9] 8 - (9 + fog o 


Design conditions, however, are usually specified in terms of gross 
weight and load factor rather than root angle of attack. The inclusion 
of these two additional independent variables requires two additional 
equations. By considering that the airplane must be in equilibrium as 
regards vertical forces and pitching moments, the two additional equa- 
tions may be written as 


2|2h| {7} + Pp - nW = 0 (17) 


for equilibrium of vertical forces and 
-2|2nx| 2) + ea [nc?| feng} - PiX + nWx, = 0 (18) 


for equilibrium of pitching moments about the pitch axis. 


Equations (16), (17), and (18) are the basic equations for a 
flexible-wing airplane. They may be solved simultaneously for {i}, 
Œp, and Pp as functions of any design values of speed, gross weight, 
and load factor. Equation (16) as written applies to symmetrical flight 
conditions, but by substitution of an antisymmetrical El matrix 


(see appendix A) it is then applicable to unsymmetrical flight condi- 
tions which are considered in greater detail under the following sec- 
tion heading. 
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Equations (16), (17), and (18) can be altered as shown in the fol- 
lowing paragraphs to include (a) the effects of forces on the fuselage, 
(b) the effects of external stores by the method derived in appendix D, 
and (c) the interference effects on the exposed-wing loading due to the 
presence of the fuselage by the method derived in appendix E. 


In connection with forces on the fuselage (item (8)), the lift and 
moment characteristics of the fuselage in the presence of the wing are 


assumed to be known. The lift and pitching moment of the fuselage may 
then be written with small error as 


ee E b (cra), | (19) 


Mp = qSc Pen + (oy), (20) 


in which fuselage up loads and nose-up moments are considered positive. 
This lift and moment may be appropriately included in equations (16), 
(17), and (18) to get the following more complete set of equations 

(see fig. 2): 


For the wing load distribution, 
O 
EE Ela m 
for the summtion of vertical forces, 


2|2n] t + as (Cur) es + Pp = OW - 38r. (22) 
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and for the summation of pitching moments, 
al |) + Eb (0), a = (Cp), 0 + Pm = 


nW = + 2q EN + (meg - Xp) (23) 


The twist term Po] appearing in equation (21) has been included 


at this point to provide for the possibility that tail loads may enter 
the wing at some point along the span, as for a tail-boom type of con- 


figuration, for instance. This Pip Of contribution is otherwise con- 


sidered to be zero. A method for handling the tail-boom type of airplane 
as well as the case of the tailless airplane is described in appendix F. 


In considering the effects of the external stores (alteration (b)), 
as in the case of the fuselage, the Lift and moment characteristics of 
the stores in the presence of the wing are assumed to be known so that 
the lift and moment can be given by expressions similar to equations (19) 
and (20) for the fuselage. 


The lift-distribution equation for the store case derived in appen- 
dix D, with the Pp term included as before for generality, is 


El = a% m) E - E {7} + Pp{o} - for} = 
{og} + a8 8} (2%) 


in which the matrices n and {s} relate an a, type of twist of 


the elastic wing to the store lift and moment. This twist is handled in 


two parts so that the part that varies with the span loading {2} may 
be introduced on the left-hand side of the equation in a manner parallel 


to that for the [s2] matrix. 
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The vertical-force and pitching-moment equations (22) and (23) are 
modified to include store 1ift and moment in a manner similar to that 
in which the fuselage lift and moment were previously included. In the 
notation introduced in &ppendix D these equations for one or more stores 
become 


2 |en| + 298 n c Pal Ex E ul , 
as (Cop), + Pp = nW - QSCi, - 298 $ Lo, (25) 


for sumation of vertical forces and 


2|2hx | + 298 y (Ci), [En || - > (en). Pal E 2l hj + 
aSx (Cr, ) or - 950 (Cp) or + Prop = 2a enc? | (emo) - 


MS tng, DE-UL 3S (Cp - CI.) + nWx, - Iy9 (26) 


for summation of pitching moments. In these equations the subseript n 
is the store index which, for any particular configuration, has as many : 


values as there are stores on the semispan. The matrix notation [En | 
which gives the angle of attack of the store is defined in appendix D. 


The third and final alteration (item (c)) to be considered in con- 
nection with equations (16), (17), and (18) is the change in span load 
distribution of the exposed wing due to the presence of the fuselage. 

This interference effect is derived in appendix E. The over-all fuselage- 
interference effect is considered to be composed of two parts. The first 
is that due to the image vortex system within the fuselage which is 
required to satisfy the condition of zero velocity normal to and at the 
fuselage surface. This condition is satisfied by adding a correcting 


JH 


x eer -—- 
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EA matrix (see appendix E) to the [81] downwash matrix. The second 


part of the effect is the increment in vertical velocities over the 
exposed wing due to the presence of the fuselage at an angle of attack. 


This effect is calculated as an interference twist of the les] type 


and is expressed as a function of fuselage angle of attack times the 


O 
"overvelocity" matrix [So]: (See appendix E.) The lift distribution 


(eq. (21)) when altered to include these fuselage-interference effects 
becomes 


Eis-spo-psdoe: 
Cs) ] «[- Y (27) 


O 
where the elements of the [So | matrix give the increments in vertical 


velocities along the span. 


The calculation of these fuselage effects would not be required if 
appropriate data were available from wind-tunnel tests of a scaled model 
of the subject airplane. A method of determining these and other aero- 
dynamic twists as well as the applicable values of section lift-curve 
slope from appropriate wind-tunnel data is given in appendix G. The 
method utilizes equation (12) to obtain aerodynamic coefficients which 
are free of model wing flexibility effects and which are therefore appli- 
cable to the full-scale airplane having a wing flexibility different from 
that of the model. 


Unsymmetrical flight conditions.- In addition to the symmetrical 
flight conditions already outlined, a number of unsymmetrical flight 


conditions are usually investigated in structural design. Among the 
conditions which may readily be investigated by the methods of this 
report are those which arise through the use of roll-producing devices 
such as ailerons or spoilers. The load distributions on an elastic wing 
associated with roll-control deflections may be thought of as the suma- 
tion of distributions from the following specific loadings: 


~ 


re ii e cmm eee ym p e n e RS A 
t - P aoe 4 aA" PLE EN LI am x 
* a M ’ A + " =. Pp ES = 
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(1) The symmetrical or pre-roll loading which existed prior to the 
deflection of the roll-producing device. 


(2) The incremental loading resulting from the deflection of the 
roll-producing device. 


(5) The incremental loading associated with a constant rolling 
velocity with no roll-control deflection, which is usually described 
as the damping-in-roll loading. 


() The incremental loading caused by the rolling angular accelera- 
tion. This loading results from the structural twists {ag} of the 


elastic wing when the rolling acceleration p acts on the wing mass 
distribution (see appendix C). Note that the resulting aerodynamic 
rolling moment will generally be in the same direction as the applied 
rolling acceleration. 


Three flight rolling conditions will be used to outline the pro- 
cedure for determining unsymmetrical loadings on the Wing. These 
rolling conditions are: 


(a) Steady roll at some specified value of wing-tip helix angle pb/2V 
with no rolling acceleration 


(b) Roll initiation resulting from the instantaneous deflection of 
the roll-control device to the angle required to obtain the specified 
value of pb/2V but with no rolling velocity 


(c) Roll terminstion, that is, control surfaces deflected in 
opposition to the steady rolling velocity  pb/2V 


Steady roll: In a steady-roll condition the span load distribution 
for the elastic wing is given by the summation of the first three loadings 
enumerated. 


The distribution obtained for the first, or pre-roll, loading. (is 
described in detail in the section Meine trical Plight conditions." 
Equations (16), (17), and (18) or their appropriate equivalents are 
used together with the symmetrical [51] matrix of appendix A. 


The distribution of loading and its associated rolling moment Ml 
or Xs, resulting from a unit antisymmetrical deflection Da] or 65g) 
of the control device are obtained from equation (16), the antisym- 
metrical matrix [31] of appendix A, and the appropriate {ag} matrices 
which give the aerodynamic and structural twists resulting from a unit 


deflection of the roll-control device. These fag} matrices can be 


written in terms of control-surface deflection by the method described 
in appendix F or by means of data from wind-tunnel tests, if available. 


———— F —- — - a a et RR Ree are y n eg n Aa DA A SS RT a i m y p ———Mm———— À————À — —— - 
N - rS a N " 2 - 
» 2 af - r man - 
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Similarly, the distribution of loading and the damping moment Mo 
associated with a unit value of the wing-tip helix angle (pb/2V), are 


obtained from equation (16) and the antisymmetrical [s1] matrix. The 


{ag} values in this case vary linearly and antisymmetrically across the 


span from (pb/2V), at one tip to -(pb/2V), at the other. 


These unit load distributions associated with ô or ô and 


al 81 
(pb/2V), must then be scaled up or down according to the amount of con- 
trol deflection Sg or 6, required to give the desired value of pb /2V. 


The deflection required is obtained from the equation of equilibrium of 
the &irplane in roll as 


Sa] ev (pb/2V)4 
or 
Oo M 


. NER. UE 
os, &V (pb/ev), 


where the rolling moment produced by the control deflection balances 
the rolling moment due to damping in roll in the steady specified rolling 
condition. 


After the unit load distributions have been scaled in the manner 
just described, they may be added to the pre-roll loading to obtain the 
load distribution for the specified steady-roll condition. 


In this outline the assumption of equal and opposite roll-control 
deflections is made. If, as is more generally the case, unequal deflec- 
tions of the ailerons or spoilers are involved, the span loading must be 
determined in a slightly different way. To illustrate the procedure, a 
spoiler deflection on only one wing may be considered to be equivalent 
to a symmetrical and an antisymmetrical deflection with an amplitude 


- eee = - — -— + — er e ER rr n a e A rr M 99 ii 
x j id " H * a LI " . a n uS Es ` ee ^ i e 
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equal to half of the spoiler deflection 6,/2. The antisymmetrical 


deflection results in a gain in lift on one wing which is exactly bal- 
anced by & loss in lift on the other wing so that & rolling moment is 
produced without a change in the over-all wing lift. The symmetrical 
deflection, however, results in a change in total wing lift, and hence 
load factor, with no rolling moment. In order to compensate for this 
change, both the vertical-force and the pitching-moment balance equa- 
tions must be introduced so that the wing load distribution associated 
with the change in a, required to compensate for the change in tail 


load can be determined. Adding the symmetrical and unsymmetrical loading 
distributions will give the proper lift distribution for the wing with 
a single spoiler deflection Og. 


Differentially operated ailerons might be considered in & manner 
similar to that outlined for the spoiler. A further extension to include 
the combined deflection of both ailerons and spoilers can also be made 
in the same way by making use of the proper spoiler-to-aileron gearing 
ratio. 


Roll initiation: In the roll-initiation condition where no rolling 
velocity is assumed to exist, all the listed loadings occur except that 
due to damping in roll. 


Since the control deflection (8g or 65g) will already be known 


from the steady-roll condition, tbe problem is to determine the initial 
rolling angular &cceleration p for instantaneous control deflection. 
The procedure involved is first to find the wing spanwise airload dis- 
tribution and its rolling moment ES due to a unit rolling angular 


acceleration P4- The values of ME. depend on the wing mass and 


stiffness distributions as well as wing aerodynamics. The values are 
obtained from equation (16), or its equivalent, for antisymmetric flight. 
With the value of Mp. known, the desired angular acceleration p is 


then found by solving the following equation of motion for the airplane 
in roll: 


rl m (298) 
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or 
Ms 
5g T 
p = —— (29b) 
g 
386.4 Py 


where the moments are given in inch-pounds. The value of p obtained 
from equation (29) is then used to scale the loadings previously found 
for the unit rolling acceleration to the correct value. 


Roll termination: In the roll-termination condition, the airplane 
is assumed to be rolling with a wing-tip helix angle pb/2V and the 
roll-control device is moved abruptly in a direction such as to reduce 
the rolling velocity to zero. As in the roll-initiation condition, the 
desired airplane rolling acceleration is obtained from the equation 
expressing the equilibrium of the airload and inertia-load rolling 
moments, which for aileron control is 


The airloads on the wing are those caused by the pre-roll condition 
plus the airloads from aileron deflection, damping in roll, and rolling 
acceleration. The inertia loads are those arising from the pre-roll 
condition plus the effects of the rolling acceleration jp. 


DISCUSSION 


The method outlined in this report not only includes several previ- 
ously omitted items which are of practical interest in the design of a 
wing for aeroelastic effects but also is sufficiently extensive in scope 
that almost any type of airplane configuration may be considered. Because 
of its length, however, the method is better adapted to the determination 


of loads on & specific airplane rather than to preliminary design studies 
of several configurations. 
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Matrix formulation of the problem has particular merit for such a 
general treatment since discontinuities in angles -or masses due to either 
special aerodynamic or structural features can readily be included. It 
also permits future improvements to be made to the details of the separate 
appendixes without essentially modifying the method which has been outlined. 
In fact, throughout the paper it is possible that the engineer would modify 
the method to suit his own needs and draw from the appendixes whatever mate- 
rial would be required to investigate the problem at hand. In such a modi- 
fication he would of necessity consider the relative merits of ease in com- 
putation against the accuracy both of the method and of the data available. 


For these reasons, only a few general guides which might be considered 
for successful application of the method are given. 


For mg equal to 2s, equation (12) will give essentially the 
same results as those given by the Weissinger I-method of reference 2 
which is valid for wings of arbitrary plan form and having flat-plate, 
circular-arc, or parabolically cambered airfoil sections (refs. 5, 4, 
and 5). The method would be expected to give the most accurate results 
when applicable values of my are known, such as those obtainable by 
the method of appendix G for instance, since in general the fuselage, 
flaps, and external stores will affect the applicable values of mo. 
in fact, equivalent values of any of the aerodynamic parameters as 
obtained from experimental data by the method of appendix G are pre- 
ferable to purely theoretical values and may easily be incorporated. 


The treatment of compressibility effects used in this report, 
Wherein each wing section is permitted to have its own compressibility 
correction, differs from the Prandtl-Glauert method in that the wing 
plan form is not distorted; instead, the angles of attack are altered 
as indicated by equation (11). The treatment adopted has the merit of 
considerable saving in time for equal or better accuracy since only 
one 1] matrix is required for all Mach numbers. The methods of 
obtai compressible values of my are described in appendixes A 
and G. 


With regard to the number and selection of the horseshoe vortices 
to be used to represent span loading, it is suggested that the horse- 
shoes be chosen narrower over that portion of the span where large gradi- 
ents in loading are expected, that is, near the ends of control surfaces, 
near large changes in sweep, and at the wing tips. At least two vortices 
should be used with each control surface and a minimum of seven per semi- 
span is suggested for a ‘clean’ wing. 


With regard to the structural parameters EI and GJ required, it 
may be stated that equivalent values which include the stiffnesses con- 
tributed by the leading- and trailing-edge structure should be used in 
preference to the usually conservative values employed in the structural 
analysis of the wing for shear and bending stresses. 
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In the design of a simple wing-fuselage combination without external 
stores or nacelles, it is necessary to study and apply the results of 


only appendixes A, B, C, and E which describe the |s|, [So | ; and foz} 


matrices as well as the fuselage-interference effects. Effects of store 
and nacelle loads and. moments on the wing load distribution are covered 

in appenđix D. Appendix F outlines the modifications required to adapt 
the method to the determination of wing loads on flexible tailless and 
tail-boom airplane configurations. Appendix G indicates a procedure for 
obtaining equivalent values of section lift-curve slopes, effective com- 
pressible section sweep angles, and interference twists from wind-tunnel 
tests of models which may not be scaled correctly for flexibility. Appen- 
dix H deals with the determination of the divergence dynamic pressure of a 
swept-wing airplane with a large external store. The problem of divergence 
normally does not occur with a swept wing except that the attachment of 

a large external store may cause it to diverge. The determination of 

the divergence dynamic pressure is the only case in this paper which 
requires iterative procedures. 


Boeing Airplane Company, 
Seattle, Wash., July 8, 1953. 
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APPENDIX A 


AERODYNAMIC FUNDAMENTALS 


The [51] Matrix 


The purpose of this appendix is to explain, in simplified form, the 
aerodynamics involved in and the steps necessary for the computation of 


the downwash matrix EN and to develop the correction for compressi- 
bility used in this report. 


The lift or circulation distribution which varies along the span 
of a wing can be visualized as resulting from a system of horseshoe 
vortices, each of which is of constant strength. Such a system of 
horseshoes is illustrated in the following sketches, in which double 
arrows are used:to indicate that the sense of circulation around each 
line-vortex segment is given by the right-hand rule: 


la Ip Ie la Le lf 


Actual airload 
curve 


43 —— = Approximation to the 


actual loading as 
given by horseshoe 
vortices 
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It is obvious from the sketches that the shape of the actual load 
distribution may be approximated to any practical degree of accuracy by 
a suitable change in the number of horseshoes, each of constant strength. 
The point of importance is that the net strength of the trailing vortex 
at any point on the span of the wing is numerically equal to the rate 
of change of strength of the bound vortex in the spanwise direction. 

The strength of the trailing vortices would therefore be greater for 
those portions of the wing span over which the more rapid changes in 
the spanwise airload distribution occur. 


Results of theoretical investigations have shown that little loss 
in accuracy with respect to the spanwise airload distribution will be 
entailed if: 


(1) The total strength of the chordwise system of bound vortices 
is concentrated in one bound vortex located at the local streamwise 
quarter-chord point. 


(2) The downwash angle at each vortex station across the span of 
the wing, at a point one-half of the local streamwise chord downstream 
of the bound vortex, is equal to the geometric angle of attack.  Herein- 
after, this point is referred to as the downwash control point D. 


The downwash angle at any such control point D is therefore the 
total induced downwash velocity at that point, normal to the plane of 
the wing and caused by the complete system of bound and trailing vortices, 
divided by the flight velocity of the wing. 


It should be mentioned that the condition described in paragraph (2) 
is true as written only for airfoils having a two-dimensional lift-curve 
slope equal to 2r. As is discussed in a subsequent section of this 
appendix, the condition described in paragraph (2) is modified to the 
following form when the section two-dimensional lift-curve slope is 
different from 21: 


O, Es os 


The downwash velocity at a point due to a single horseshoe which 
is composed of a bound vortex and its associated pair of trailing vortices 
is known to be proportional to the strength of the circulation of that 
horseshoe and, therefore, proportional to the running lift on that por- 
tion of wing span represented by the bound vortex of that horseshoe. The 
downwash angle at any one downwash control point thus is the sum of the 
incremental downwash angles due to each of the horseshoes in the system 
of horseshoes which represent the wing and its Lift distribution. 
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Within the limitations which &ccrue from the use of & limited number 
of horseshoe vortices to represent & wing, the elements with which to 
solve the following problem are now available: 


Given the geometry of a wing plan form, the angle-of-attack varia- 
tion, and the section two-dimensional lift-curve-slope variation across 
the span, determine the spanwise airload distribution. The unknowns 
are, of course, the values of the rumning lift at each of several points 
on the span. A necessary condition to the determination of these unknowns 
is that as many independent equations be available as there are unknowns. 
This condition can be fulfilled if the angle of attack is known at each 
of the wing stations for which the loading is to be determined. 


It is obvious that, if the strength of each bound vortex represents 
the average airload over its own portion of the wing span, good accuracy 
will be obtained if the values of the running load, as determined from 
the solution of the simultaneous equations, are considered to be valid 
at the midspan point of each bound vortex. The shape and distribution 
of the continuously varying airload curve is then obtained by fairing 
a curve through all of the points thus obtained, with the restriction 
that the loading must drop to zero at the wing tip. 


The method for determining the EM matrix'is now illustrated for 


a typical wing shown in the following sketch, which includes a system 
of horseshoe vortices and associated downwash control points: 


V 


1 
Np Typical 


_-Leading edge 


—- Locus of c/4 points 


a 
po 
~ 


Trailing 
edge 


These aa 
Trailing vortices "Locus of downwash control 
extend to infinity points at 3c/4 
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In actual practice it will usually be necessary to use a larger number 
of horseshoes to represent the wing. Four per semispan are used in 

the following presentation for convenience only. The following typical 
information is available at each spanwise station: 


The section lift-curve slopes (in two-dimensional flow) at sta- 
tions m - Tn; Mb. . . are mg = M; m5, . ... . The angles of attack 


of the section zero-lift line at stations 1 = Tq» Nor - + - are 
Lp = 04, Q5, . e- 
Since & linear relationship exists between the strength I'3 of a 


particular horseshoe vortex 'j and the downwash velocity ME ata 


particular point i on the wing plan form due to that horseshoe vortex, 
the following general equation can be written: 


w.. -K,.D (A2) 


where K is a constant. A particular horseshoe vortex, such as I, 


then causes the following values of downwash velocity at control points 1 
to 8: 


Wi] = Eu Woy = Koli 
Woy = Eo Py Wer = Seals 

(A3) 
Vig a a Ya “at 


Similar relationships exist between To, lx, ... I4 and the control 
points 1 to 8, that is, 


ME = als Vos = K l5 e o e è o Wgz = 8515 (ad) 


26 NACA TN 3030 


If the horseshoes are symmetrical with respect to the center line 
of the wing plan form, the following values of K are equal for the 
typical wing sketched: 


Ki; 7 Kgg — Kip = Kg; — Kig-Kgg Xi = Kgg 
Koi = K7g — Kop = Kyy Koz = K76 Ko, = K75 
Kai = K68 Kao = Égq Kzz = Kgg Kay = Kes 
Sei 558, “at AGS 596 M= 
(A5) 
Kor = Kg Bsa = Eq  Ks3=K6 —— Ko, = Kys 
Koi = Kag Kéz = Kaz Kg = K365 Ky, = K35 
Ki 7 Kog = Kya = Koy Kp = Kog K7} = Kas 
Re =%ig Ben = %7  Eg=kEg Ee = 5 


Further, for symmetrical wing plan forms at zero yaw it is always 
possible to obtain a spanwise airload distribution which is unsymmetrical 
with respect to the plan-form center line as the sum of two airload dis- 
tributions, one of which is symmetrical and the other of which is antisym- 
metrical with respect to the plan-form center line; for example, 
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At each spanwise station the S component is that due to the symmetrical 
distribution of load, the A component is that due to the antisymmetrical 
distribution, and the U component represents the algebraic sum of the 

S and A components, that is, the unsymmetrical distribution of load 
across the span. This division offers a considerable reduction in the 
amount of work required in that, for either symmetrical or unsymmetrical 
flight conditions, airload distributions need be determined on only one- 
half of the wing, provided, of course, that proper account is taken of the 
sign of the circulations existing over the other half of the wing. 


For & symmetrical distribution of &irload over the span 


(A6) 
lo = I4 Py = I5 
and for an antisymmetrical distribution of airload 
D4 = Tg X = T6 
(A7) 


The total downwash velocity at any control point is the sum of the 
downwash-velocity contributions at that point that are induced by each 
of the horseshoe vortices in the system that represents the wing; that 
is, 


TW 


13 + "uy + Wig + Met up V 


T 18 


(48) 


Wo = Wa] + Wop + Woz + Woy + Wos Wa + Wog + Wog 
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If the equation for the downwash velocity at control point 1 is 
expanded as an example, the result is 


Wi = Ki + Kioo + Ks t K yy t Kis + Kj, eG t 


+ Egg (49) 


Alternatively, by use of the relations in equation (A5), equation (A9) 
can be written as 


Wy = Ki + KL, t Ms 4 + KE, T Kats + Eg + 
Kgol4 4 Kgil'g (A10) 


In case a symmetrical airload condition is being investigated, sub- 
stituting equations (A6) into equation (A9) results in 


We = Py (Ky + K 8) + is + Kj) T PX) + K1g) t 


Ty, Ky, + Kis) (A11) 


For an antisymmetrical airload distribution, substituting equations (AT) 
into equation (A9) results in 


Ty (Ki - Kis) (A12) 
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For the typical horseshoe system assumed, the complete series of 
equations relating downwash velocities to the circulations are: 
For a symmetrical airload distribution, 
Wis = TK + Kg) + Too + Kig) rs (5 + Kug) + 
Py (Ky, + K15) 
Was = Py (Kp, + Kog) + To (Kos + Koq) + T'3(Kog + Kog) + 


Ty (Koy + Kos) 


(A13) 
Wa, = Py (Ks) + Kg) + Po(Xzo Kap) + Ps (Kzs + Ka) + 
Py, (Ka, + 1,5) 
Wig = Taya + Kya) + Pollo + Ej) + T5 (5 + Kye) + 
ry (Ky + Kys) 
and for an antisymmetrical airload distribution, 
nag r (K E Ky g) + (Ka - E 7) + Pz (Ky, - Kg) + 
Py (Ky, - Ky5) 
Wa, = T1(Kor - Kog) + To(Koo - Ko7) + "(Kos - Kog) + 
Py (Koi - Kos) (A1) 


Wy = Ti - BQ) + lo (Ko - Koz) + "s(&55 - 159) + 
Pus - K55) 

Wip = (By. - Eig) + Po(Kyo - Kyr) + P3(%,3 - Kyo) + 
Py (Kyy - Kis) 
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From equation (Al), which expresses the relationship that must exist 
between the downwash angle w/V at each control point, the wing angle 
of attack ap, and the section lift-curve slope m, for the wing station 


at the control point, the following series of equations result, where, 
typically, Wy represents either Yig or En 


Vo 2n V Ont 
(A15). 
ta "a Hh TT 
vo 2x vo x 
If 
kj] = An(Ky, +t Kg) 
kjo = 4a (Kyp + Kj-) 
yy = Wt (Kyy + Kys) 
or 


where the upper sign is used for symmetrical airloads, the lower sign 
is used for antisymmetrical airioads, and the subscripts L and R in 
equation (A16) mean left and right wing, respectively, then substitution 
of equation (A16) into equations (Alt) gives 

hw, = KQQDI4 + KQ5T5 + ki5Iz + Ky 1, 

to = koj] + koola + ko3P3 + Kon lh 


haw = Xy 1 + Kolo + IL + Ky} Ih 


(A17) 
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The relation between the running load 1 and the circulation is 


(A18) 
hy = eva, 


and the equation relating dynamic pressure to mass density and true 
airspeed is 


q = pv" /2 (A19) 


Equations (A15), (A17), (A18), and (A19) may be combined to give the 
following final system of linear equations: 


ky + Kolo + Ejzdg + Kyyty = ham 
Koy ty + Kool + Koss + Koy ty = Hamoa 
(A20) 
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Equations (A20) can now be written in matrix form. The equation 
for each station can be written as 


y 
2 Ky yt, = hama 
j=l 


n 
ECUEEOT 
j=l 
(A21) 


h 
NEUE 
$-1 


h 
> kj ¿24 = hamo, 
$i 


Equations (A21) and therefore equations (A20) can be expressed in the 
general form 


- 
> E = hqm,a, (i =1,2,3,...n) 


In matrix notation equation (A22) becomes 


630) - pegiz aen 
where 
E11 Kip kis E 
Koy Koo Koz Koy " 
e 
[Si] Kz] Ezo kzz Ez) (a2 ) 


Kyl Kyo Ky Ey) 


———— ——— ——— — —- -——— — - - 
oe 
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by 
41) = “2 (A25) 
bs 
y 
and 
uo, 
O Lama 
DNE " leu (A26) 
ham) on 


The development of the steps necessary to compute the Kij elements 
in the EI or "downwash" matrix follow. As a first step the relation 


which exists between the strength of & segment of & straight-line vortex 
and the induced velocity at a nearby point should be found. If the 

strength of the vortex (whose sense is given by the right-hand rule for 
moments) is T, the velocity induced at the point P can be written as 


. T(eos a - cos B) 


where a and f are the angles between the direction of the vortex seg- 
ment and lines joining the ends of the segment to the point as shown in 
the sketch 

End 


D, inZ/sec T 
e 
a 
view 


R, in. 
“oF >, —8 => Wp, in./sec 


Plan view 


Equation (A27) is used in the subsequent derivation of the ET matrix. 
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A plan view of the geometry of a typical horseshoe vortex on the 
left-hand wing is given in the following sketch, in which distances and 
angles are considered positive.as indicated and the sense of the circu- 
lation of the elements of the horseshoe vortex is given by the right- 


hand moment rule: 
E 


YD 


| -r - — "V 


D | 
l $ 
Sx 
V 
D-- 
Sy . ` 
. l These dimensions ore 
those of the actual 


| Vortices of strength wing plan form 


I extend to infinity 


@ Dimension . 
reference point 


The points V and D are a typical horseshoe reference point and a 
typical control point, respectively. 


The incremental downwash velocities induced by a single horseshoe 
vortex, if downwash velocities are considered as positive, are: 
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(1) For the left-hand trailing vortex the relations 
R = -(s, - h) 


a= 07 cosa = l 


B 


270° - 60 cos B = -sin 9 


are substituted into the general relation 


_ T(cos a - cos B) 


W. (A28) 
a lc 
to obtain the incremental downwash velocity 
(l + sin 0) 
"^T Wee, - i) ed 
Ey - 5) 
(2) From the right-hand trailing vortex where 
R = Sy + h 
a = 90° - g cos a = sin Ø 
B = 180° cos B = -1 
the incremental downwash velocity is 
l(sin $ + 1) 
gute. (430) 


hx (8; + h) 
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(3) From the bound vortex where 


R = 8, 
a = 180° - 9 cos a = -cos O0 
B = 180% -ø cos B = -cos Ø 


the incremental downwash velocity is 


Tr(cos Ø - cos 0) 
liste. 


(A51) 


The total downwash velocity at a typical control point due to a 


complete single horseshoe is then 


pf Asem’, rm mot 


By - h Sy + h By 


Substitution of the identities 


p {L+ sing nie) 
W = =] e e crm 
cos Ø cos 9 


(432) 


(A352) 


(A34) 
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Equations (A34) and (A2) are identical; therefore, 


(A55) 


If, as indicated in the sketch 


Typical horseshoe y =y 
referencè point. 


Sx 


\ . 
Typical 
control point | 


n 


the control points are assumed to be located on the left semispan of 
the wing so that 6; and gr. represent the pertinent angles for a 


horseshoe located on the left semispan and Op and fr represent the 


pertinent angles for the corresponding horseshoe on the right semispan, 
then for a typical control point 


Iia E 


cos Gr, cos 6, 
(A36) 


1 + sin fj -l + sin 0g 
CM cos Pe cos Og 
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From equation (A16), 


k = bt (Kr, £ Kg) 


so that in the [$1] or downwash matrix each element k is computed 
from the equation 


i fi + sin fy L+ sin, 1+ sin fe _1+ sin Op 


Fs a ee Py cee 
SN cos f, eos 81, cos Gp cos Op 


(A37) 


As in equation (A16), the upper sign is to be used for symmetrical span- 
wise airload distributions, and the lower sign is to be used for anti- 
symmetrical spanwise airload distributions. Note that in equation (A3) 


ae DY 

(Xp - Xy) + (p - Xy + b)? 
cos Ø = d. Ji un 

(% y) + (vp - xy + 1)? 

DN 

sin 0 = 

(% - x) + oyy 
cos 0 = zi NE ils 


)H 
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The ET matrix therefore is computed from the matrix equation 


- (A38) 


ES 


where the upper sign is used for symmetric flight conditions and the 
lower sign is used for antisymmetrical conditions. 


1 E sin fr, l+ sinó, 1+ sin Øp l + sin Og 
MI + A s 


cos PL cos 0, T COB Ép * eos 9 


Since the EN matrix is used in equation (A22), which is 


EJH- Bex] 


the elements of the. [S1] matrix are seen to be influence coefficients 


relating the incremental downwash angle at each control point to the 
intensity of the running lift over each increment of the semispan of the 
wing. In general, all the elements in the principal diagonal of the 


EN matrix will always be positive and those elements not in the 


principal diagonal.will always be negative because the velocities were 
considered as positive downward and wash velocities from a horseshoe 
vortex are downward only in the region behind the bound vortex and 
between the trailing vortices of that horseshoe. 


Compressibility Corrections 


The method by which compressibility effects are handled in this 
report is based on simple sweep theory. This theory is presented in 
references 6 to.8 and substantiated in references 9 to 12. Summarized 
briefly, infinite-aspect-ratio sweep theory indicates that compressi- 
bility effects are functions of the effective Mach number Moa, which 
in this case is the Mach number of the stream velocity component that 
is normal to the leading edge, so that the lift-curve slope in com- 
pressible flow is given by 


In 


lig eerie (459) 


V1 - M?cos?Aw 
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where 
m section lift-curve slope at M= 0 
M stream Mach number 


effective sweep angle for compressibility effects or yaw 
angle of infinite wing, radians 


The same relationship exists for the effects of Mach number on 
section pitching-moment coefficient 


ET: SA (Aho) 


Cn, = 
yı - M°cos“Ay 


Compressibility effects on a finite-aspect-ratio swept wing can be 
handled in a similar fashion. For the finite-aspect-ratio swept wing, 
however, recognition should be given to the fact that compressibility 
effects will vary across the span of the wing. In general, smaller 
values of the effective sweep angie are indicated for the wing root and 
tip sections than for the midsemispan region. Even in the midsemispan 
region the effective sweep angles for compressibility effects are not 
functions of the wing plan-form geometry alone; the spanwise variations 
of camber, thickness ratio, chordwise thickness distribution, and angle 
of attack are likewise involved in the determination of the value for 
the local effective sweep angles. 


Consider equation (7) for the downwash angle induced at a distance 
rearward of the lifting line which was derived from two-dimensional 
considerations 


Wr 
V 


d 


i 
2p ^f 


This equation expresses the relation between the downwash angle at the 
three-quarter-chord point (i.e., at r - c/2) and the geometric angle 
of attack ap when the section lift-curve slope is equal to its theo- 
retical value 21. 


Several approaches can now be made to the compressibility problem 
regarding the values of lift-curve slope, the angles of attack, and the 
location of control points to be utilized for equating downwash angle 
to geometric angle of attack. In the method used in this paper the 
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distance r is always equal to one-half of the local streamwise chord, 
so that in matrix notation equation (7) reduces to (see eq. (9)) 


t) 3c/4 H t] 


The applicable values of mg to be used in this equation are those 

at the effective Mach number Me = M cos Ay, and thus each wing station 
is permitted to have its own compressibility corrections. In the absence 
of test data reduced according to the method of appendix G, a value of 
m, can be calculated from equation (A29) for a value of Am equal to 


the geometric sweep angle of the plan-form streamwise quarter-chord line. 
Although, in general, this procedure will yield only an approximation 

to the correct theoretical value, this value will usually be sufficiently 
accurate for preliminary design purposes when A/cos A is large (where 
A is the aspect ratio) or M cos A is small, or both. 


This "modified angle of attack” method is characteristically dif- 
ferent from the "plan form distortion methods of handling compressi- 
bility effects (refs. 13 and 14); the former distorts the angle at the 
original three-quarter-chord point, whereas the latter stretches the 
half-chord. distance. 


An important advantage in the compressibility-correction method 
presented in this paper is that only one [8 | matrix need be computed, 
whereas a different (8 | matrix for each Mach number is required in 
the plan-form-distortion method. 


The following developments show the equivalence of these two methods 
of accounting for compressibility effects. The first comparison is for 
the case of an infinite wing in yaw. 


If AM is the sweep or yaw angle for this infinite wing, then 
according to the plan-form-distortion procedure (indicated by the sub- 
script pd) the equivalent wing in incompressible flow (M = 0) is to 
be at a yaw angle such that 


= _ tan Ay 


A - 32 


Further, the lift-curve slope of this equivalent wing in incompressible 


flow, multiplied by the factor 1/y1 - M^, is equal to the lift-curve 
slope of the original wing in compressible flow. 


tan 
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The lift equation for the equivalent Mach number method of this 
paper (indicated by the subscript Me) is 


Im, = qca E 
VA = M°cos“Ay 


and that according to the plan-form-distortion theory is 


2x cos Ay 
CL = 


fi - we 


In order for the two lifts to be equal, the following equality mist be 
true for any value of M: 


Ind = qc 


cos Ay _ CoS AM 


‘al - M ! yl - MÉcos?AM 


It can be shown that such is the case by the use of elementary trigo- 
nometric identities. 


The next comparison is for the case with zero sweepback. The plan- 


form-distortion method stretches the half-chord distance as a function 
of Mach number to 


l 


Vi - ME 


which when substituted for r in equation (3) gives 


! 
c = 
2 


MIA 


NACA TN 3030 43 


Equating the downwash angle at the stretched control point to the 
section angle of attack gives 


Bu ue (A1) 
V icy 


The compressible lift-curve slope (ref. 14) is 


mo = —£— (a42) 


Ji - e 


Equations (A41) and (Ale) may be substituted into equation (2) to give 


em rVi- M? py?c 
yl = Me xcV 2 


or 
l = pvr 
For the compressibility-effects method of this paper the downwash 


angle at the three-quarter-chord point equated to the angle of attack 
corrected for effects of section lift-curve slope is (see eq. (8)) 


Y O A: 
(t). yh xecV en of 


from which 


ap = E (A15) 
mocV 


When equation (A43) is substituted into equation (2) the same final 
answer is obtained: 


t= pvr 


The identities of compressibility corrections in the case of either 
the unyawed (Ay = 0) infinite wing or the yawed (Ay # 0) infinite wing 
substantiates, at least qualitatively, the use of the compressibility 
correction methods as presented in this paper. 


ee i ey A AAA A AAA A ee — ni ig ee AAA 
, * pe c" ou A a 
» * A = 4, s+ 


mp m e tee Pe am ana A AN ee ee a 
dJ * 
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APPENDIX B 
THE ELASTICITY MATRICES [S5] AND [S,'] 


The purpose of this appendix is to explain, in brief form, the 
structural theory involved in and the steps necessary for the computa- 


tion of the elasticity matrices [S| and [s37]. 


Development of the ES Matrix 


In the development of the downwash matrix El in appendix A, a 


continuously varying spanwise airload distribution was replaced by a 
series of constant-intensity running loads, each of which covers an 
increment of the wing semispan. This concept of an equivalent system 
of loads is also used in this appendix. 


Consider the geometry of the structural skeleton of the wing to be 
as represented in figure 5, in which double arrows indicate that the 


right-hand rule of moments applies. The notation in figure 3 is defined 
as follows: 


L4, Lo, -- - GL, total lift of increment of wing having span of eh, , 
eho, . . . 2h,, respectively, numbered inboard 
from the left wing tip, lb (see eq. (B1)) 


e streamwise distance from horseshoe reference point 
at a wing station to elastic axis at the same 
station, positive when elastic-axis point is to 
rear of horseshoe reference point, in. 


f streamwise distance from elastic-axis reference 
point at a station to elastic-axis reference 
point at: the next inboard station, positive when 
inboard elastic-axis point is forward of outboard 
elastic-axis- point 


d5; e.. lateral distance between elastic-axis point at a 
station and elastic-axis point at the next station 
inboard 


17? 


A,» Mo» y see A, local sweepback angle of elastic axis 
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Mx, > M. 


MY: MY? ; 


PEE 


15505 es 


The general 


Ti 
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rolling moment et elastic-axis point around longi- 
tudinal axis through local elastic-axis reference 
point due to total lift of all the vortices out- 
board of this point, positive when it raises left 
wing tip, in-1b 


pitching moment at elastic-axis point around lateral 
axis through local elastic-axis reference point 
due to total lift of all the vortices outboard of 
this point, positive when nose up, in-lb 


beam bending moment at elastic-axis point about an 
axis perpendicular to local elastic axis, positive 
when it puts compression in wing upper surface, 
in-lb (see eq. (B2)) 


torsional moment around elastic axis at local 
elastic-axis point, positive when it is in direc- 
tion of leading edge up, in-1b (see eq. (B3)) 


form for the wing lift is 


Ly = 2bghy (m1) 


where ly is the intensity of the running lift at station n measured 
in pounds per inch and 2h, is the span of the horseshoe vortex at 
station M3 the general form for the bending moment is 


My = Mg cos An - My sin Ay (B2) 


and the general form for the torsional moment is 


- monem A A HA AAA AAA m a e mm T e a —À —- —À 
LI . ` * re n a . 


Ta = My. cos A, + My sin A, (B3) 
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At station 1, on the center line of the horseshoe vortex nearest the 
left wing tip (see figs. 3 and 4), the following equations apply: 


My, = EL TL (Bi) 
1 2 2 n 
I1 / hy tan Ay (95) 
Wy = BRL 2 
Which, when substituted into equations (B2) and (B5), yield 
hi e- sin Ay 
M = L,(————_- - a (B6) 
4 cos A 2 
e, cos A. 
1 1 


At station 2, on the center line of the second horseshoe vortex 
inboard. of the left wing tip, 


Lobo 
Mz, = I4, + ko (B8) 


Llo[ = bp ten Ap - 
erate +e ET) o 
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and equations (B2) and (B4) become 


Mo = Lydy cos Ap - I4 (ei - fi)sin Ap + eg. - A (B10) 


cos Ap , 2 
T5 = L,d, sin A, + I4 (e, - f) coB A, + V [87 ^e) (B11) 
At station 3, 
h 
My, = Ly (A, + dp) + Loa, + za (B12) 


and, similarly, | 
Ms = [E (8, + ap) + La, cos A - mo -fj-£j)s 


h ez sin 
Lo (ee x fa) sin Az + of < m E uw (211) 


oe | (87 m do) * Lodo] sin Ag + Ha - t] - £5) + 


ez m 


" ta - tie Az + Ly ————- (2315) 
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Equations similar to these can be developed for M and T for 
the remaining stations on the left semispan of the wing. 


In algebraic form, the equations for M, and T) are 


n=l n-l n-l n-l 
DEED ELEME RO 
k-l l-k k=1 t=k 


by ey sin Ay 
: cos A, y 2 (B16) 
n-l - n-l n-1. 
T, = sin Ap > (Ix d, + cos An > Lar (Ex - > f + 
k=l t=k k=1 t=k 
dod (B17) 
2 


Notice that the first two terms in each of these equations are equal 
to zero when n - 1, that is, for the tip station. 


In matrix notation, equations (B16) and (B17) become 


{mp = [= A][r,] - es à] e (B18) 


{ql - ps A] [ro] + qe A] B (B19) 
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where 


M 
Mp 
(1) = 4M, 
Mn 
ni 
To 
{t} = 43 
T, 
| In 
Lo 
i = 35 
La 
cos Aj O O O 
O COS Ao O 0 
0 cos Az O 
O 0 cos Ay 
sin A1 O Q 0) 
O sin Ap 0 O 
O sin Az, O 


T T ` - 
, "m E] a xd + = 
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(B20) 


(B21) 


(B22) 


(B23) 


(B2)) 


b “a ` Li ` “ly 
# bi t - a a - 
- ^ . >» = 
js » 
Da A 
d d 


E E eo - =- 4 - -r >” - 
^, i : E - e QC CU EE P ME gue, A à IN 
U vs ` - re » e E e r * + Ll - = 
T , ` r 2M. .* . PEE . 
TE E ` -7R ee ur e - & . 
- * ~ " E " ` 
ne aa a AA A A A tl al lr en — gl. © a 


F] = 


0 O 
e - £4 O 
e, -f, - £ eo - f5 
e1 - fi =to- f3 eo - fp - f3 


di + dp 


dj + da + à; 


do + ds 


do + d; + dy 


0 
0 
0 
E de 
0 
0 
0 
by _ en tan Ay 
4 cose A, 2 
dy 


O 0 
O 0) 
0 O 
O 0 
ey ~ fh, 0 

0 

0 

0 

0 

EN RENS each 

l costs e 


(B25) 


(B26) 


OS 


OCOC NL VOVN 


= - à = * - ~ e > icd . : : : g a E 
" - « a z A 
E . » « = be ` y AN "m 
` . P axi ub wu P * Pe > 2 A 4 i . 
T * E Ri < n . ‘ 
lll — — ——— M— te ee E 
m — H»'"——— si i a m e a i i II i i Ed A tla Pe sccm vical AE uc qe ECRIRE TE PE: S a an — We aai — À oo pa rar 


O O O . 
2 tan Aj 
dy 2 0 0 ; i 
2 tan Ao 
EE c 209 0 
di + do do 8 tan As 0 
T5| = e. - | (B27) 

[Fe] dj + da + dz do + dz da 2 tan Ay, 0 ° . 
dy + dp + dz + d, da+dz+ dy dg + d, dy, Y C 

. an 


Equations (B18) and (B19) provide the means for computing the bending and torsional moments 
along the span of the wing due to the loads L as given by equation (Bl). 


Each of the loads 2 or L, however, is affected by the variation in angle of attack from 
root to tip. This problem is handled by computing the streamwise angle-of-attack change Oe 


due to the structural deformations caused by the series of loads (ij Since a streamwise angle- 
of-attack change is required, assume that a unit positive pitching moment my (nose-up moment) 


is applied in turn at each of the points 1, 2, ...n on the elastic axis. The unit pitching 
moment is in the plane in which a, is to be measured and its direction agrees with the posi- 


tive sense of Qs. ' 


OSOS NL VOVN 


TS 
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The streamwise angle-of-attack change can then be obtained by 
applying the general relation (see ref. 15) 


o 
mM ds tT ds 
ABg 7 i +f (B28) 
O EL 0 GJ 
where 
sy desired streamwise angle-of-attack change due to values of 
bending moments M and torsional moments T along elastic 
axis of wing resulting from the series of loads inj 
m beam bending moment per unit pitching moment applied st sta- 
tion at which ag, is to be determined 
t torsional moment around elastic axis per unit pitching 
moment 
ds increment along elastic axis 
BI effective beam bending stiffness around the axis of the bending 
moments M and m 
GJ effective torsional stiffness around the axis of the torsional 


moments T and t 


The stations on the wing for which the angle-of-attack changes «a, 
are to be computed are those on the center line of each horseshoe vortex. 


It is assumed that sufficient accuracy in the results will be 
obtained by using the values of bending moment M and torsional 
moment T obtained by means of equations (B18) and (B19) and that these 
values of M and T may be considered to be constant over the incre- 
ments of the wing span corresponding to the span (2h) of each horseshoe 
vortex. 


Values of EI and GJ are also assumed to be constant over each 
such increment in wing span, and the values to be used &re the effective 
values which correspond to the wing station y at the midpoint of each 
horseshoe vortex. 


The general relation (eq. (B28)) can be used to express a, in 
matrix form by letting 


- ————— — —— A - CoC M eee am Rt mur" que ciet? cm nm: sud" dun LL O ge C d: QD. EM coe er — m de —À— Du 
7 - E , - ^ : —á = = = : a "- = - T = fl ee " - —— —-— y ~a 
"o -". * ` : * - f ` t» ae a et, e . ^. hl y, e 4 E " ^ = 
Soe T * € ~ = ` ` - . n e E: 2 we "X * + à = z " 
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An sweepback angle of elastic axis at station 1, radians 
MA beam bending moment at station Tj, due to loads {i} on wing 
Ta torsional moment around elastic axis at station 1, due to 


loads {i} on wing 


D. j beam bending moment at station i per unit pitching moment 
applied at station j 


torsional moment around elastic axis at station i per unit 


: pitching moment applied at station j 
2hn 
(às), = ——- 
COS An 
eh, horseshoe span at station Th 
AN angle-of-attack change at station Ty due to all the 


loads (1) on wing 


Then, for station 1, 


may Pha Man 


¿y (EI), cos A (EI), cos Ao o (EI), cos Aq 
t414174h 2t54T ol 
Lo E id ME IL Ec 2% —— mmm (B29) 
(GJ); cos Ay (GJ)o cos Ap (GJ), cos Ay 


for station 2, 


am Mh Mealy tain 


Qs 7 


(EI), cos A, (EI)5 cos Ap E (EI), cos A, 

2t oTihy topToho thot hy 
HA ee Se abere (B30) 
(GJ); cos A, (GI)5 cos Ay (GJ), cos A, 


and so forth. 
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From the geometry of the problem, when point i is at or inboard 
of point j, the bending moment due to a unit pitching moment applied at 
point j is 


my; = -sin Ay (B318) 


when point i is outboard of point j, the bending moment is 


mi, = 0 (B31b) 


when point i is at or inboard of point j, the torsional moment due to 
a unit pitching moment applied at point j is 


t,, = cos Ay (B32a) 


and when point i is outboard of point j, the torsional moment is 


a 4 = 0 (852b) 


J 


and the complete series of equations for «az, of which equations (B29) 
and (B30) are representative examples, can be written in matrix form as 


ae cde o 
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Typical elements in [m] and [t] are 


-sin A;/2 -sin Ap -sin Aj -sin Ay 


O -gin Np [2 -gin Az -sin Ay 
O O -8in A /? -sin A, > : 
=p. . (253) 
O O -sin Ay, /? ` " 


cos A; /2 cos Ao cos Az cos Ay 
O cos Ap fe cos As co8 Ay 
0 O cos Az /? cos Ay 


O O O cos Ay/2 


(B35) 


[+] = 


: Substituting equations (B1), (B18), and (B19) into equation (B33) 
gives 


fas} - Be = 3 = ap] - e tasi: 


E an cos a xs A][ro] + [os A] s] exp) (B36) 


ee ee AA LD ÁÓ v e Home dpa E ru m Pe itn ct 
2 » a. , » 
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if 


[sg = (62 E Jl es JE - pes g| + 


cos =] B ES Ala) + Io > Afi Ei] — Gen 


then 


fe - 226) e 


Equation (B47) represents the most general form for the [So] matrix, 


and each element &4 4 of this matrix represents the angle-of-attack change 
in radians at station i due to the structural deflection of the wing 


caused by a unit loading at station j. In effect, the [So] matrix is 


an array of influence coefficients, and the elements of this matrix may 
be computed according to equation (B37), or, when an actual wing is avail- 
able, they may be obtained by load-deflection tests of that wing. 


Development of the Auxiliary Elasticity Matrix ES 


The EJ matrix as defined in the previous section considers air- 


loads outboard of each specified station, as defined by a streamwise cut 
through the wing box. Inasmuch as ribs are conventionally installed in 
a plane essentially normal to the elastic axis, the spar~box loads are 
more correctly determined by considering the airloads outboard of a plane 


LI 
Ld y E ve 7 - v we on e mn a a a gp gi i al A EE A ae A cuan — £r rh Re. A mm amet me —— eee 2 
- weet g mi m ae + - "A ` * * ” 
g H T > . Ju M T & ola! > ^ € as * * 
‘ * - ` , - = ^ r h 
e <= m ` P cm E . A E P , b - e- 
a > ‘ ^. 009 4a e xn ^ ros = . a 1 
- + a “eto at s "E d PERLE a 
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normal to the elastic axis. These loadings are represented in the fol- 
lowing sketch: 


Trailing- 
edge 


Nn 


The corrective loading is then the loading obtained by subtracting 
the loading of (A) from that of . The triangular areas (C) and O) 
represent the corrective loading. It will be noticed that the loading 
over the area (C) is the negative of the loading of the corresponding 
area in À and the loading over (D) is the positive of the corresponding 
area in . ln order to obtain these corrective loadings, it will be 
necessary to assume a distribution of pressure over the corrective areas. 


The equations are now developed that represent the shear, wing 
moment, and wing torsion due to the corrective loadings, and the 


EA matrix is modified to include these effects. A plan view of the 
left wing is shown in figure 5, where 


E local angle between elastic axis and lateral axis equal 
to A, the sweep angle of the elastic axis 

L local angle between leading edge and lateral axis 

T local angie between trailing edge and lateral axis 


The correction terms for moment, shear, and torsion due to wing 
airload in the EX matrix are 


os — en on -——— - 
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OAM 3 aoe correction to rolling moment at elastic-axis point 
i around longitudinal axis through local elastic-axis 
E AMy reference point, positive when it raises left wing 
tip, in-1b 


Mr , AM , correction to pitching moment at elastic-axis point 
i 2 around lateral axis through local elastic-axis 


9 a Ay reference point, positive when nose up, in-lb 
AM), AM, correction to wing bending moment at elastic-axis point 
e. o. AM, about an axis perpendicular to local elastic axis, 


positive when it puts compression in wing upper 
surface, in-lb (see eq. (B39)) 


AT, Ato, correction to torsional moment around elastic axis at 
local elastic-axis point, positive when it is in direc- 


B tion of leading edge up, in-1b (see eq. (BhO)) 
A54 , So, correction to local shear, positive up, lb 
o LÀ e AS 
n 


The wing bending moment M, is 


My = My, cos E, - My sin E, (B39) 
and the wing torsional moment a is 
T = My cos E, + My sin E, (BLO) 


Assumed pressure distribution.- The distribution of pressure p 


over the corrective areas is assumed to be given by the two-dimensional 
equation for & thin flat plate 


p = haga (B41) 


and modified to account for the actual value of local running lift. In 
the modified form, the distribution of pressure (measured in pounds per 
square inch) along a streamwise chord becomes 


(Bh2) 


nen a — I p co e rm c pU Roc mof ` 
* a vus ES " = p X M. nto 
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Where x/c is the fraction of the streamwise chord to the rear of the 


leading edge. The ratio 2/c is assumed to be constant so that the lines 


of constant pressure radiate from the origin O in figure 5. In accord- 
ance with the assumptions previously given concerning the corrective 
loadings, the pressures are negative (downward) over area FDG and posi- 
tive (upward) over area ABG. 


À polar coordin&te system (0,0) is selected, and the shear, moment, 
and torsion due to the pressures on area OFB will be subtracted from 
the shear, moment, and torsion due to the pressures on area ODA. Lines 


of 0 = Constant are lines of constant pressure and p will be a func- 
tion of 8 only. 


From the geometry of figure 5, p may be expressed as a function 
of 0 by writing 


c + OM tan T - x = OM ten(L - 0) 


Hence 
x= c+ [tan T - tan(L - 0)| (B43) 
and 
c 4 OM tan T = OM tan L 
or 
= c 
DUAE. (Bh) 


tan L - tan T 


Combining equations (B43) and (Bhh) yields 


(B45) 


x= cjl + 


tan T ~ tan(L - 6) 
tan L - tan T 


If equation (B45) is substituted into equation (Bh2), 
- - T 
p=k tan(L ~ 6) - tan (Bh6) 
tan L - tan(L - 6) . 


x22 (B47) 


where 
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Rolling-moment correction at station n, AMX,,:7 The rolling-moment 
correction at station n, AM? is obtained by adding the rolling 
moment AMy, due to area OFB to the rolling moment AMen due to 
area ODA (fig. 5). 


The rolling moment AMX a due to area OFB may be found by writing 


the equation for the elementary moment about FM for any value of p 


between O and that at the line FM as: 


dAMy = pom - p cos(L - SIE dp dO 
The equation of the line FM is 
p cos(L - 0) = OM 
or the value of 'p at the line FM is 


OM 
cos(L - 0) 


Over area OFB the moment is negative; hence, 


- OM 
— 6r cos L-0) 
AMy, = -KOM tan(L - 9) - tan T ag p dp + 
A E tan L - tan(L - 0) O 


OM 
L-T cos(L-8) 
x tan(L - 6) - tan T cos(L - 9)ao ¿246 
O tan L - tan(L - 0) 0 


Since O0 varies between O and L - T, the integrals can be evaluated 
by letting 


tan L - tan(L - 9) 2 sin2u 
tan L - tan T 


| 


(B48) 


ten(L - 0) - tan T | | 2j 
tan L - tan(L - 0) 


20 SS D eee pcm um des Ce n ————— MÀ — M — — —— - a tl 
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lt follows that 


a(L - 0) = -2(tan L - tan T)cos^(L - 0)sin U cos U aU 


Consequently, 


E x/2 
AM. = -KOM (tan L - tan T) cos“U dU + 
O 


" f? 
E OM? (tan L - tan T) E cos^U dU 


z 2 KOW (tan Lo tan?) 


102 
POP: MT (B49) 
6(tan L - tan T)? 


Ihe rolling moment AN due to area ODA may be derived in a 
similar manner to get 


L-T 


nom = E cos ORT 
Mf, = KOM ten(L - 0) - ten T tan T a9 duds 
B ten L - tan(L - 0) 9) ^ 
R 
COSl L-0-E 
ten(L - 0) - tan T y cos(L - 0)à0 p^àp 
tan L - a - 
O 
(B50) 
where 
R = OM cos(E - T) + e(1 - €)sin E (B51) 


cos T 


me ee OCU —— — 
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After integration this equation simplified to 


xKOMRe(tan L - tan T) 
AK, = M - 
„(tan L + cot E)^sin?E|i - o? 


xKR^(tan L - tan T)(h - 308) 


2h(tan L + cot £)?stz^g (1 - q2)? 


(B52) 


where 


2 _ tan L- tan T 
tan L + cot E. 


(o Sa <1) 


Equations (Bh9) and (B52) may now be combined to give the rolling- 
moment correction AMy,, at station n: 


AMy = ony, + 5 
xk 6cR?(1 - Q°)(tan L + cot E)sin E - RÓ(% - 30%) (tan L - tan T) 
24 (a - 02) (tan L + cot E)?sin?E 
2c? 

eine | ees RO (B53) 

(tan L - tan T)? 
If U defined in equations (Bh8) is written as 

Uel¢ cos(E - T)cos L (854) 


sin E sin(L - T) 
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and if 
V = tan L + cot E (B55) 


then R (see eq. (B51)) becomes 


R = c(U - e)sin E (B56) 


6(u - el (1 - q2) - (U - e)0(h - 3Q2)Q2 2 


= —- a ____—_ B57) 
JE Pa eya An 


For the special case of the untapered wing Where the leading edge is 
parallel to the trailing edge or T = L, 


2 
_ tee sin°E cos^L 


e? = he + 1) (B58) 
E cos?(E - L) 


Pitching-moment correction at station n, AMy .- The incremental 


pitching-moment corrections due to the areas OFB and ODA may be 
derived and combined in a manner similar to that for the rolling-moment 
corrections and.the equation for AMr, becomes 


;c? | CU - 9)? |(} - 302) tan L-eN|«ea- q2) (U- e) (eQ2v - tan L) 


n 


E ! ve(i - 92)3/? 


(6e -1)VQ°-2 tan L 


(859) 
y2q} 


For the special case of the untapered wing where the leading edge is 
parallel to the trailing edge or T = L, 


e 
te 2 sin B cos L sin E cos L tan L 
My = - —— (8e° - 4 pepe aee remi B60) 
n 16 DT cos(E - L) 1 cos(E - L) i 


mt a a m i, nome ce tiem n 
* 
r 
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Shear correction at station n, ASpn.- The elementary shear is 


dASn = pp dp de (B61) 


From this equation the incremental shears over the areas OFB and 
ODA may also be derived and combined in a manner similar to that for 
the rolling moment to give 


Sy = =|. - — (B62) 


For the special case of the untapered wing where the leading edge is 
parallel to the trailing edge or T= L, 


ic sin E cos L j J 
E 4 cos(E - L) "m (B63) 


Mođification of En matrix.- The Ea matrix is modified by 


using the corrective rolling and pitching equations. Using the equations 
= - in 
AM, = AMy cos E, - AM sin E, (B6) 


AT, = AMy sin E, + Aly cos E, (B65) 


and letting (from eqs. (B57) and (B59)) 


-AMy,,/2 (B66) 


T 
| 


J = CAM; (B67) 
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results in 


AM, = (-p cos E, + J sin Ep)? 


AT, = -(p sin E, + J cos Ep)? 


or, in matrix form (with A, = Ej), 
{ost} = E E fein A] EO, (B68) 
Lam) = - E Apu] + [os 6] {i} (B69) 


Equation (B18) thus becomes 


fu} - Ee 3 (61123 - 6| - E» 9/0183 - a] f) Gro 


and equation (B19) becomes 


fx) - E a Eae] - 03) + ee does -i DEC 


The corrected form of the [S2] matrix thus becomes 


1 - EST nas e) e an - + 


pea E gleac- 6] + e a pate - si 


(B72) | 
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H 


where, as in equation (B48), 


(as) = [82423 (875) 


It should be noted that, because of the assumed pressure distribu- 
tion given by equation (Bhl), this correction does not correctly handle 
conditions involving deflected flaps and roll-control devices. 
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APPENDIX C 


COMPUTATION OF ls 3 MATRICES 


This appendix outlines the method by which the component parts of 
the fag) matrices may be determined. ‘The (2%) matrices comprise two 


classes of twist: (I) those which would be present even if the wing were 
rigid and (II) those due to inertia effects, thrust or drag, and section 
pitching moments on the flexible wing. A given airplane design condition 
may obviously require any combination of the twists listed under classes I 


and II, that is, 
{ag} = tes]. 4 led. (C1) 


For simplicity, the sources of fag) may be summarized as follows: 


Class I - aerodynamic twists (i.e., zero-lift-line shifts due to 
effects other than wing structural deflections) 


(a) Built-in twist due to camber or construction, or both. 


(b) Interference twist due to fuselage, external stores, and 
so forth. 


(c) Twist due to flap deflection. 

(d) Twist due to aileron deflection. 

(e) Twist due to spoiler deflection. 

(f) Apparent twist due to airplane rolling velocity. 


(g) Apparent twist due to airplane pitching velocity. Angles of 
attack due to airplane pitching velocity should be measured at 3c/h. 


Class II - structural twists due to wing deflections caused by the 
following inertia and aerodynamic loadings which are independent of the 
wing lift distribution: 


(a) Vertical acceleration acting upon dry-wing dead weight, 
wing internal-fuel dead weight, and external-stores dead weight. 


(b) Effect of airplane rolling acceleration upon dry-wing dead 
weight, wing internal-fuel dead weight, and external-stores dead 
weight. i 
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(c) Effect of airplane pitching acceleration upon dry-wing 
dead weight, wing internal-fuel dead weight, and external-stores 
dead weight. 


(d) External-stores thrust or drag. 


(e) Section pitching-moment coefficient with flaps, ailerons, 
spoilers, and so forth in neutral position. 


(f) Incremental section pitching-moment coefficient due to flap 
deflection, aileron deflection, or spoiler deflection, or in various 
combinations. 


Of the class I twists, the aerodynamic built-in twist is known from 
the wing geometry and the characteristics of the profiles used, and the 
interference twist may be obtained from actual tests of the complete air- 
plane configuration, either model or full-scale, by using the procedure 
of appendix G or from other calculation methods. Twists due to control 
deflection are determinable from the same type of tests as were used to 
evaluate the interference twist. Apparent twists due to airplane rolling 
and pitching velocities are completely determined when the airplane flight 
condition to be investigated is specified. 


The type of twists due to the effects of wing deflections arising 
from loads which are independent of wing angle of attack, such as those 
listed under class II, may be computed with the aid of equation (B33) 


Since, in this case, [as] z [ss] 


e BEA OBB e 


where jd and [rl are the wing bending moments and the torsions &long 
the wing elastic axis due to the loadings of class II or to any combina- 
tion of them. 


Although the Tg) twists do not all explicitly require the cal- 
culation of loadings, they do influence the equilibrium airload and must 
be accounted for in determining the net wing loads from equations (21), 
(22), and (23). 


NACA TN 3030 69 


APPENDIX D 
DERIVATION OF EXTERNAL-STORE MATRICES 


The term external stores as used in this report is intended to apply 
to such items as nacelles, external fuel tanks, bombs, rockets, and 
similar items commonly attached to the wings of airplanes. 


In this appendix the lift coefficients and pitching-moment coeffi- 
cients for each external store in the presence of the rest of the air- 
plane configuration are assumed to be known and to vary linearly with 
angle of attack. On this basis a set of linear equations is developed 
in which the airloads on each external store are accounted for in com- 
puting the deflections of and airloads on the elastic wing or airplane. 
The coefficients used are based on the airplane wing area and the wing 
mean aerodynamic chord. There is a specific lateral reference axis for 
each external-store pitching-moment coefficient, and the angle of attack 
of the external store is taken to be the angle of attack of that wing 
station within the span of whose horseshoe vortex the station of the 
external store lies. Special symbols used in this appendix are as 
follows: 


Cr? C; lift and pitching-moment coefficients, respectively, of an 


external store measured at its aerodynamic center (see 
eqs. (D1)) 


- 
Che , Cg lift and pitching-moment coefficients, respectively, of an 
O O external store at zero angle of attack measured at its 


aerodynamic center 


Imo» Mao lift and pitching moment, respectively, of an external store 
measured at its aerodynamic center (see eqs. (D2)) 


My , My ; rolling moment &t elastic-axis point around longitudinal 
€ etl axis through local elastic-axis point due to loads on & 
: My n particular store outboard of this point, positive when 

it tends to raise left wing tip, in-1b 


; pitching moment at elastic-axis reference point around 
s eri lateral axis through this reference point and due to 

: Mya loads on & particular store outboard of this point, 
positive when it tends to raise nose, in-lb 


Additional symbols required are shown in the following sketch: 


———— ee ee ee Mn 
E m Pe m t a ttp D SS | LLL m 
i M Aem - ^ a ^ o s A EE ZEE NER Oe D pn ng 

E ` Ld . +” ., i , : 
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Forward 

,Aerodynamic cenfer 
i of external store Inboard 
i Ye 
uL 

+ M Me +2 
| o zao Te+2 

Tt ec? 
Xe E+l 


+ My 
| 
|. he ` 
Consider one external store whose characteristics are 
(D1) 
Cm = Cmpa T (Cm) tg 
So that 
La = 38 [Org + (Prog | 
(D2) 


Mao = aS [ens + (Ca) | 


LOH 
NACA TN 3030 71 


The lift and moment on the external store Lig and Mao are located 


with respect to point e by the dimensions y¿ and xg as shown in 
the sketch. 


At the first station inboard of the external store €, 


Mx. = Love 
| My 7 Mao + LgoXe 
Me = Mp cos Ac - My, sin Ag 


= Lgo(ye cos Ae - Xe Sin Ac) - Mpo sin Ac (D3) 


Te = My. Sin Ac + My. cos Ae 
=Imo(ye sin Ae + xe cos Ac) + Mpo cos Ae (DÀ) 


At the next inboard station € + 1, 


My rl = Imo (Ye + de) 


eun = Mao + Lgo(xe - fe) 
Me+1 = LEO Gre + dejcos Mey, ~ (xe - fc)sin hess | - Mgo sin Ae) (D5) 


Te+1 7 Eo m + dc)sin Ag; + (Xe - F¿)cos Aq, | + Mag COS Agcy, (D6) 


. ==. 
a 
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At station € + 2 


J 


My erp z Lro (Ye + de + de+1) 


My. = Mro + Lgo(Xe - fe - ferl) 


Mero = Igo Gre + de + de,1)e08 Agro - (xc - fc - fc,1)sin Merz | 7 


Mgo Sin Aero (D7) 


Te+o = Igo (ye + de + de+1)sin Ae+2 + (Xe - fe - fey])cos here | + 


Mao cos Nero (D8) 


At station n, the most inboard station considered, 


n-l 
E PED 


Since 


My = Mx, COS An - My sin Ay 


In = Mx, Sin An + My, cos An 
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then, 
n-l / n-l 

Mn = Imo |cos AnlYe + 3 - gin a(t = > ta) - Mao sin An 
In-€ m=€ 


(D9) 


and 
n-1 n-l 
Tn = Lgo [sin taye + > dm} + cos bale - > 3 + Mag cos Ay 


m=e€ 


i 


(D10) 


In matrix notation, the following equations, typical for each 
external store on the semispan, are equivalent to the algebraic equa- 
tions (D9) and (D10): 


fa) = SERA 8) - eda] (d foun - d] Guo uo 
(sj 


where, if the value at the top of the column is for the tip and the 

other values are for the successive stations inboard to the root and 
ce is the station within whose horseshoe span the external store is 

located, 


li 


(Es A] R} - [eos A] (9); Lgo + [cos A] T Meo (D12) 


O 
O 
O 
€ 


fr} = d (D13) 
Ye +t de 
Ye + de + de+1 
Ye t de + detl + dejo 


- 


weer —- ———À a cay mc nn " am Meme y TOTEM onim cce i Ri a D MM Pg — UI i M P MÀ X Ó — —— —— — M — MÀ M — — 
` ? e « - 
» 1 
B 
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ooo 


tvi = Xe - fe (D14) 
Xe - fe - fal 
Xe - fe - fetl - fe+2 


(D15) 


fs} 


Pp|Éoocd 


{Ms} ^ Meu (D16) 


{te} = 4 Te (D17) 


From equation (B33) 


E ls a os 


— — — — 


. 
- = "ye z s 
Toe o er ey 2 " - 
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where the external-store bending moment and torsion are given by equa- 
tions (D11) and (D12), respectively. 


Substitution of equations (D2) into equations (D11) and (D12) leads 
to: 


ts} 7 (ee alf} - [ein Dens, + {eo alg} I 
[sin A] e (Cr) see - [sin A is basen - 
[sin a] (tg) WBE (Cn, ) ag | (019) 


and 


fog} = de A18) + ES a] pastry, + ess i « 
Eds fi} hs (Gg nts + Eo? A] Exjastiny, + 
eos ^] fe Jase (eng) arty , (9) 
— T 
£a) = fiet AI fe} - [3 ie) (21) 
fron} = (ls a} + ES 163) (022) 
then equations (D19) and (D20) become 
{Ma} = as Teo, + fKou}(Crg) of - [Ein A] fs] Song, - 
[sin a] fig} Cag) (p23) 
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fi} = ee fena, * (rer) eae, + Eo? io], + 


[cos A] Bg) (Cm) 2 T" 
Now let 
fen] 7 t ex jCts - [sin A] fig}, (D25) 
Kou} = om) (Gr) <7 [sin A] Ug) e (Cg), (D26) 
and 
fa] = orp CE + Cos a] En pECmg (D27) 
{Kor} = [Xov (eng), + [eos a] Exjo(Cmg)a, (D28) 
then 


{Ma} = 48 m * (edens | (D29) 
f = 96 {fas} + bera, y (D30) 


Equation (D18) relating {os pp with {Ma} and (Tg ) can now be 


written as 


{cs | = q5 1] + fen (D31) 


= -—,— = * 
ae * a D = — Pe os Oe nir e 0Á— —— i 
a * . AAA 
4 , -— on . 
. 


- æ nd ——— ———M— 
. .. . 
" . è EI a E 
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where 


fa} = ug x| A PE Bedae ia} 
al) En ied 


The equations just derived are typical for a single external store 
on the semispan. Each additional external store on the semispan requires 
a similar set of equations. 


(D32) 


The total angle-of-attack change due to the wing airloads shear and 
the airload on the external store is 


where the subscript W denotes wing and E denotes external store. 


Here there is an Josh term on the right-hand side for each exter- 
nal store; that is, 


[ssl = {ost + [s]. + fag 2, Pk» (D34) 


Equations (12), (13), anā (B38) for {i}, fas}, and fag} „p respectively, 


repeated for convenience are 


AE e 


f) = {ox} * (08) + 9] 
[S2] Ln (95), 
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Including the effects of external stores on {ag} according to 
equations (D34) yields 


E Se] CF = food + feet Rute, teu s oom 
From equation (D31) 


1%] Ep pahe + ij ES: 


(D36) 
{estn = aaas + {ahoen 
However, ap, can be found from 
arp = |Enjler) (D37) 
where a typical value for [En | ls given by 
je] o o 20000 c (a 


The row matrix | En | has as many columns as there are horseshoe 


vortices on the semispan and has a value of unity entered in the column 
which corresponds to that horseshoe within whose span the external store 
is located. (In the example implied in eq. (D38) there are ten stations 
and the external store n is within the span of the fourth horseshoe 
from the left wing tip.) 


Combining equations (12) and (D37) results in 


«s "sf ea id 


LIH 
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Combining equations (D35), (D36), (D39), and (B38) yields 
PIO = fb + foc + Bs > 
qs fe, + Eeh ET ass [s1] { a 


O 


fed. oed] 


4- 


+ 
E 


[4 y [e]. [En | EN z is (D41) 


[9-2 E] E xa 


If equations (D41) and (Dh2) are substituted into equation (D40), 
the following matrix equation is obtained for an elastic wing having a 
number of external stores on the semispan: 


[: - q8 nj e EN - i] {2} = {ap} + deg + asf}  (Dk5) 


This equation is similar in importance and usefulness to equa- 
tion (16) except that the aeroelastic effects of the airloads upon the 
external stores are now included. It reduces to equation (16) if all 


the [4] and {B} matrices are set equal to zero, since [1] = [1.6]. 


- —— n_e e o ———— e —— — o —— We — "o - ——— MM e MM, cnc mio Án 4M MÀ —Ó— MM m ———— —— — «— ——— A A ———————— — - — = 
- = ` z 
; M e’ ` " , * - a . = 
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APPENDIX E 
WING-FUSELAGE INTERFERENCE 


In this appendix an approximate method of calculating the effect 
of a fuselage on the spanwise airload distribution on the wing is devel- 
oped. The primary sources of this interference are considered (ref. 16) 
to be as follows: 


(1) The effects of the image vortex system inside the fuselage on 
the downwash angles at the various control points over the wing semispan. 
This image vortex system induces a flow which is a first approximation 
to that necessary to satisfy the condition that there be zero velocity 
normal to the fuselage. 


(2) The effects of the vertical ' overvelocities" resulting from the 
fuselage angle of attack. These velocity increments affect the local 
angle of attack at the various control points on the wing. 


This method of analysis can be used when applicable data of the type 
deseribed in appendix G are not available. 


Item (1) may be considered as the effects of wing airloads upon 
themselves due to the presence of the Puselage; whereas item (2) accounts 
for the effects resulting from the fuselage having an angle'of attack. 


In the development that follows the fuselage is assumed to be of 
circular cross section, of constant diameter, and infinitely long. 


For the image vortex system mentioned in item (1)., the individual 
images of the wing trailing vortices can be shown to be located on a 
straight line joining the axis of the fuselage with the axis of the par- 
ticular wing trailing vortex at a distance from the fuselage center line 
such that 


Ro = E (E1) 


where a is the fuselage radius and R is the distance from the fuse- 
lage axis to the trailing vortex (see fig. 6). 


In similar fashion the bound vortex is assumed to have an image 
within the fuselage cross section. The image of the bound vortex is 
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assumed to lie on a straight line joining the forward ends of the image 
trailing vortices. The forward ends of these image trailing vortices 

are assumed to lie in the same transverse plane (perpendicular to the 
airplane center line) that contains the particular bound vortex being 
represented. Figure 6(a) shows a transverse section for the high-midwing 
configuration and figure 6(b) shows a plan view for the midwing 
configuration. 


This system of real and image vortices is only an approximation. 
It satisfies the condition of no flow across the fuselage boundary only 
at the transverse plane containing a bound vortex and its image and the 
transverse plane infinitely far behind the wing. In general, some flow 
will occur across the fuselage boundaries, and hence, to some degree at 
least, the corrective downwash velocities induced at the wing control 
points by the image vortices will be in error. 


Within these limitations it is therefore considered that wing- 
fuselage interference effects upon the wing spanwise airload distribu- 
tion can be obtained by: 

(1) Adding to the elements of the [sa | matrix for the exposed wing 
a corrective matrix EN whose elements represent the vertical downwash 


velocities induced at the various wing horseshoe control points by the 
image vortex system within the fuselage boundaries. 


Each image vortex within the fuselage is of the same strength and 
sense of circulation as its real counterpart and there are as many image 
horseshoe vortices within the fuselage as are used to represent the wing. 


(2) Adding to the twist of the wing (e a a correction due 


to the vertical overvelocity field arising from the effects of fuselage 
angle of attack. The correction applied to fag} is proportional to 
the quantity (a, - a4) where a, is the angle of attack of the wing 
root-section zero-lift line and Os is the incidence angle between the 
wing root-section zero-lift line and the fuselage center line. 


The above corrections are developed as follows. 


n ———————————————— a ee A n M M e 
k : >on gt c +. Lc 3 E camera as 
; 


ta ^ £ i 
= T Bis p Ne. ~ 
‘4 t7 . 
-- 
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The elements of the corrective matrix EN may be derived by con- 
sidering the geometry of the following sketch: 


Z (up) 


H 


f e | 
, Eo 
/ Te a 


(Left), 
Y 
-S;2' 
"Sy 
Eve. Ogie is a horizontal plane 


The horseshoe vortex shown in this sketch is considered to be one 
of the image vortices, the location of which is given by equation (E1). 
The strength of circulation r of this image vortex is equal to that 
of its real counterpart. The plane of the horseshoe wabo is at an 
angle y to the horizontal plane Ogie and the bound vortex ab lies 
in the plane Omhg. The point k may be considered to be one of the 
downwash control points as described in appendix A. 


Y 


a! r 
>. OO 
i X (rearward) 


00 


The total vertical downwash wy at the point k due to this image 
horseshoe vortex is the sum of the contributions due to each segment. By 
means of equation (A27) the contribution of the trailing vortex ao is 


7 X ' - H cos y 
Wo = El + $ e : ) (E2) 


- * ` 
Y ENNIO TO 
o RA No 
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that due to the trailing vortex bw is 


EPS e) 
and that due to the bound vortex ab is 
"Hes AE (Ek) 
eb lx D D 
where 
A= ck 
B = dk (E5) 
D = fk 


From the geometry of the sketch the distances A, B, and D are found 
from 


AS = (sy' + H cos y)? + (sz + H sin y)? 

Be = (a - H cos y)* + (sz + H sin 7)? (E6) 
2 2 1 2 

D” = 8, +t (Sy Sin y - 8, cos y) 


end the angles are obtained from 


6' = tan`? Ik. tan - D 

Pa By COB Y t 8, sin y- H 
g' = tant fE = tan” B 

Pb By. cos y + 8,' sin y + H 

(E7) 

A= ten | Sk. tan”? a 

ad s 
u= ten SE = tan ^ 4. 


€———————————— A i rR e TS NE SPE re ep em A cae Se E I A RN TE c CC C a A ea a E 
` ~ . + ^ x - = -+ 
A " . "E = = ast - - 
- ` + * 4 = > oe + , = 
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Summing the individual downwash’ contributions (eqs. (E2), (E3), and 
(Eh)) gives the total downwash at k due to one image horseshoe 


1 I T 1 
He: (Sy + H cos 7)( 1 + cos 4) (cos 9 - cos 6')s, I 


k hx "t p? 


Sy - H cos 7)(1 + cos A) 
- (58) 


B 


The downwash angles at the control points on the wing due to the 
fuselage images may then be written in matrix notation as 


fi}, "s EIE) - ¿210 > 


in which the desired correction matrix for image vortex effects is [ss | j 


the elements of: which are the locally applicable values of the function 
inside the brackets of equation (E8). The development for antisymmetrical 
conditions would be analogous to that for symmetrical conditions, the only 
change being that the sense of the vortices (and their images) on the 
left-hand side of the plane of symmetry in figure 6 would be reversed. 


The second interference effect, that due to vertical overvelocity, 
may be taken into account as shown in the following analysis: With respect 
to the fuselage center line there is a transverse component of the free- 
stream velocity V,, given by (see sketch) 


Vs n Ver = 04) 


Fuselage 


M : 
`~Root- section 
E ome zero-lift line 


Since the fuselage displaces the air in this transverse flow field, there 
are local changes in this vertical flow velocity. 


If V,' is equal to the total local vertical velocity in the pres- 
ence of the fuselage, then the overvelocity AV, is 


E OE Va (E10) 


E 
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This vertical overvelocity can be expressed as 


Ay, = E Yap - o) (511) 
y^ 


Contour maps from which the value of AV,/V, can be obtained as a 


function of vertical and lateral distances from the body center line are 
shown in figure T. 


These contours were developed from the equations for uniform flow 
past an infinitely long circular cylinder from reference 17 as follows: 


Vz' = V. e + e) (E12) 
Z 
1 
Vz 1 + aX cos 28 (E13) 
V Re 
KA 


where a is the cylinder radius, R is the radial distance from the 


axis, and € is measured from a line normal to the axis and V,. 
Thus 


Wz _ afcos 20 (B14) 
Vz R2 
From equation (Ell) the angle-of-attack increment Ax, at a control 
point due to the local overvelocity becomes 
A 


Equation (E15) in matrix notation is 


foci = [So] {ox = up (E16) 


o 
where [85] is a diagonal matrix whose elements are the locally appli- 
cable values of AV,/V,; that is, 


[o] = E (E17) 


-— — iA aa = an 
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By use of the matrices 


[52] downwash. matrix given by equation (A38) for the wing alone 
(the portion of the wing outboard of the fuselage) 
EM wing elasticity matrix given by equation (B37) 
EM image-vortex matrix given by equation (E9) 
BK 
[So | overvelocity matrix given by equation (E1T) 


and the equation 


{aa} = [5365 


the equation for wing equilibrium including fuselage-interference effects 
can be written as l 


Es] + [ss] n" - Bang fs + ag + as} + Este - zi (218) 


(62 « E] - Ed 0-0 CI] fe} - I 


(E19) 


Equation (E19) may be substituted without change for equation (21) for 
the wing load distribution. 


` + 
NAE 
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APPENDIX F 
EQUATIONS FOR TAILLESS AND TATL-BOOM ATRPLANE CONFIGURATIONS 


This appendix develops equations which, when properly inserted into 
the basic equations, allow solution for (1) a tailless airplane and 
(2) an airplane with the tail load entering the wing structure through 
tail booms. The equations are left in general form. Constant sweep 
angles and equal vortex spacing simplify the equations considerably. 


Tailiess Airplane 


For the case of a tailless airplane longitudinal balance is usually 
accomplished by the deflection of controls on the wing. These control 
deflections alter the span load distribution over the wing while con- 
tributing to the balance of the airplane, and therefore terms expressed 
as functions of the control deflection 5 will appear in the lift- 
distribution equation (eq. (21)) and the pitching-moment balance equa- 
tion (eq. (23)) in place of the Pp terms. Since the lift produced by 
the control deflection is part of the wing lift, Pp in the lift balance 


equation (eq. (22)) will be zero. 


The development of expressions in terms of the control deflection 
follows. : 


The required expression for 5 to be substituted into equation (21) 


consists of two dg components ; spa Wen is the apparent twist due to 


control deflection, and lI! which is the twist due to section pitching 
moment with control deflected. From equation (C2) 


TEE UNE 


For class II twists due to section pitching-moment coefficient 
(type (f) in appendix C), the following equations may be written for 


{Mon and {Thon in equation (F1): 
(Mio, = 79 fein A] [z5] lane? | {emo} (F2) 
{Thom q al A pg {em} (F3) 


úl 


a ee eee ae e Ry m a a oR a ~ — = — —— -- 
- n P + = 4 E 
- 4 a - 
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where 
1/2 O O O 
1 1/2 O O 
1 1 1/2 O 
I; = 
1 l 1 1/2 


For [m] and [t] in equation (F1), equations (B34) and (B35) may be 
written as 


E] = - fio] [ein a] (FL) 
[t] = [To | [cos a] (F5) 
where 
1 1 1 
1/2 1 1 
O 1/2 1 
Ig = 
O O 1/2 


Substituting equations (F2), (F3), (Fh), and (F5) into equation (Fl) yields 


fri, - BE [SEES NETS Bue» 


[js 252: e 


NACA TN 3030 


or 


Geg, de eei] + E dope co 


Equation (F6) gives the structural twist due to control deflection. If 
wind-tunnel data are not available, theoretical expressions for {em} 


in terms of the deflection 3 may be used. From reference 18, 


o 
(xl. -- psa 9o - 5 sin SAIS 
where 
00 = sie TE] 
09 = cos-i(2t - 1) 
BS 


and cp is the flap chord. Substituting these values for 06, results in 


fengys = - slog - E) - 2de(1 - £) (28 - e 


PS - £)91£8) (FT) 


The der term required can be obtained from the following expres- 
Sion for the lift produced by control deflection (ref. 18): 


{ers = 2 |x - 9, sin eol {5} 


for a two-dimensional lift-curve slope of 2x. For a two-dimensional 
lift-curve slope of mg, this equation can be written as 


(edo = Gap t+ leg 
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Substituting the values for 8o results in 


tus = EN 1 - = cos"! (et - 1) + I eve (1 - a] fs} (F8) 


O 
Dividing both sides by [mo | yields 


(e). = z È cos“l(2g - 1) + = 2de(1 - £) o} (F9) 


The total Ag term to be substituted for the Prp term in the lift- 


distribution equation is obtained by combining equations (F6), (FT), and 
(F9), to give 


A * lega, 


«698 [esa] - [Be à Eia e - 
h - = cos" (28 - e + i eG (F10) 


For the pitching-moment balance equation (23) the expression required 
to replace Pp may be obtained from the theoretical expression for flap 


pitehing-moment coefficient (eq. (FT)) as 


ae (F11) 


2hc? 
= 2 
Pp of 


‘Since in the system of equations given in matrix form (eqs. (21), 
(22), and (23)) there are two unknowns other than Pm, it would not be 
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possible to have more than one unknown in the expression to be substituted 
for Pm; in other words, the [5] matrix in equations (F10). and (F11) must. 


have a known distribution shape. let 5, be the unknown reference 


deflection; the {8} matrix may then be written as {8} - e where 
r 


(è is any chosen distribution shape. For the reference deflection at 


r 


the second station from the tip this equation becomes 


51/6, 


1.0 


15} = dry 85/er 


Oroot/9r 


For a constant deflection across the span, the distribution matrix 
becomes {i}. 


Balancing Tail Load Entering Wing Through Tail Boom 


For the case where the tail load enters the wing structure through 
& tail boom, the distribution of load over the wing will be affected 
when the tail load changes and will vary in a different way depending on 


where the tail boom enters the wing. An expression to be substituted into 
the lift-distribution equation (21) for Pp} is developed in this sec- 
tion of the appendix. 


From figure 8, the beamwise moment produced at the point P by the 
tail load entering the wing through the tail boom can be written as 


Pr Pr Pr 97 dm 
= — em sin Ap + = yp cos Ap t = 
"Pre 2 e 2 m-erl COS An 


where Ap is the sweep angle of the elastic axis at the entry section 
of the tail boom. 
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In matrix form 


(95, - 23) [2 E ¡En +2 + (er sin Ap + yp cos Ap){tp) (512) 


where [Lp] j {rr}, and [a] for six reference stations are defined as 


O O O 0 O 0 
1 O O O O 0 
1 1 0 O 0 O 
[e] - 

1 1 1 O O O 
1 1 1 1 O 0 
1 l 1 l 1 O 

oy tip 

O 

O 


1 
1 } Root 
O O O O O 
O O O O O 
o O O O O O 
[8] ~ O O O O O 
O O O detl O 
O O O O de+2 
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The equation for torque is 


[ree — T foy cos An - yp sin An) {tp} (F13) 


Substitution of equations (Fh), (F5), (F12), and (F13) into equa- 


tion (C2) gives an expression for the class II twist due to tail-boom 
entry: 


fers E 6 [zo] [e - flan J E Bm 5] i 


(er sin Ap + yp cos Ay) a - ala cos Ap - yp sin m [1] fir} 


The colum tó be substituted into equation (21) is therefore 


+29) 38 308 + 
ag tn ty + yp cos Ay + [2] [es cos ay = vp sin as (ty) (F14) 


This column will be multiplied by Pm in equation (21). Equations (22) 
and (25) remain unchanged for this configuration. 


— -—-——— — -= — —— 
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APPENDIX G 
METHOD OF REDUCING WIND-TUNNEL DATA 


Difficulty is usually experienced in applying data obtained from 
wind-tunnel tests to the design of a full-scale airplane, especially in 
those cases where of necessity the stiffness of the model differs from 
that of the full-scale airplane. The purpose of this appendix is to pre- 
sent a method of analysis by which model flexibility effects may be 
removed from the aerodynamic coefficients. 


Description of method.- The method utilizes equation (12), which is 


PJO = [band] fae} 


where 


(3 > {fo} + E) + (9j 


When equation (12) is used for computing the 1ift distribution fı} 
for a given full-scale airplane of any given flexibility, the [Sr] and 
ap} matrices used depend only on the particular configuration and the 
given flight conditions, except for certain aerodynamic-twist components 
of the fag} : 
applicable values for these Ce, components and for the [no matrix. 


These values may be determined with suitable data from wind-tunnel tests 
of a scaled model in conjunction with equation (12) as it applies to the 
model. 


matrix (see appendix C). There remain to be selected, then, 


The approach taken is usually applicable and is based on the assump- 
tion that the following data are available from wind-tunnel, tests of the 
model: (a) the spanwise variation of the section normal-force coefficient 
obtained from integration of pressure data and (b) spanwise variation of 
section chord-line angle of attack with free-stream direction obtained from 
model deflection data. These data should be available for each of several 
root-section angles of attack and at each of several Mach numbers over the 
essentially linear range of section lift coefficient. 
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A downwash matrix [sx] for the model should be computed according 


to procedures of appendix A. In general the horseshoe system should be 
identical with that which will be used in determination of the airload 
distribution on the full-scale airplane. The matrix for the model [Sul 


will then be equal to the reciprocal of the model scale factor times the 
[83 | matrix for the full-scale airplane, where the model scale factor is 
equal to the model span divided by the full-scale span. 


The final angle-of-attack matrix for the model may be written as 


Cd = [E Ed, Ghd * 69, 


= {a} cag + (A ^— (à) 


| 


where 
{a} angle of attack of root section with undisturbed stream 
(measured) 
(s built-in twist (known) 
La, 
(s) twist due to model flexibility (measured) 
M 
[ss] interference twist due to aerodynamic interference effects of 
Ip neighboring bodies (fuselage, nacelles, external stores, etc.) 
upon the wing (unknown) 
(1 sum of measured or known values 
meas 


Equation (12) for the model then becomes 


El, O - ES Bet 


Since measured values of i} are used, equation (G2) represents a 
system of independent equations (one for each reference spanwise station 
of the model) which do not require a simultaneous solution. 


me a i gm re o Ee HR mf ea - ee A AA ——— 
= 
` à , + - T 
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If {ea}, is assumed to be invariant with change in wing angle of 
b 


attack, the following matrices may be computed from data taken at two 
different root-section angles of attack: 


(ace) = (22 reas E To | ens (G3) 
[4-0 E 


The matrices (eer) and {ar} may be substituted for fa} ens 
and D in equation (G2), since the term fas}, has been eliminated, 
b 


and the applicable values of {no} can be computed from 
O 
8,0 (5) 
ej > x Hy 


The interference twist fae may now be computed from equation (G2) 
b 


by using the values of (2, computed from equation (G5) and the values of 
{a} and {a} neas measured at any root-section angle of attack, for exam- 


= ft 


The foregoing procedure indicates the simplest solution for the con- 
stants in the straight-line equations given by equations (G2). If desired, 
equations (G2) may each be solved somewhat more accurately for mo} 


and [e by & least-squares procedure which utilizes values of {7} 
b 


and {“} meas taken at several values of root-section angles of attack. 


If the interference twist a at any section is assumed to include 


& component that varies with change in some section chord-line angle of 
attack, the method is still applicable to a close approximation. In this 


case the variable interference effects will appear implicitly in {™o}- 


The approximation arises from the assumption (implied in the above 
Solution) that for a body at some spanwise station j there exists an 
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interference twist at some other station i which is proportional to the’ 
change in geometric angle at i; whereas in actuality the interference 
twist at station i is proportional to change in geometric angle at the 
body station j. 


Justification for the assumption can be shown as follows. Denote 
interference twist by qar, and the measured section chord-line angles at 


stations i and j by a; and Q4. Interference twist at station i 


due to the body at station j is written as the sum of a constant plus 
a variable twist 


(91), = (ar), + Ky jay 


where Kas is the proportionality constant. The final angle of attack 
at station i then is 


(ce), = 09 tms o). 


and equation (G2) for the spanwise station i is 


E ej. {7} cu AE ar o los 0 (ezo), ) 


However, 


therefore, 


& lel, G (ena) (e) + (er), ) 600 


The middle term is small compared to the first term and may be neglected 
since, for station i close to station j, the difference hs 7 Oy is 
negligible, and for station i far from station j the interference 


> A ar aeaa ——— are meer 
— —— ——— I NI c pc oM ber cc qM E — i a 
w . - 


a -_— n_n eee ——- HA -= —— 
e 
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effect and therefore k is negligible. If the middle term is neglected, 
equation (G6) becomes equivalent to equation (G2) and the factor 


l * k;; is apparent in the values of m, computed by equation (G5) 


even though it does not appear explicitly; that is, m, computed by 
equation (G5) is equal to the term m,(1 " ki) in equation (G6). 


Although only interference twist was considered in the foregoing 
description, the method is obviously applicable for determining other 
aerodynamic twists (see appendix C), for example, that due to flap 
deflection. 


The foregoing method of analyzing data is relatively simple and 
straightforward and has the following advantages over other currently 
available methods: 


(1) It provides a means by which data obtained from wind-tunnel tests 
of a properly instrumented elastic model wing can be evaluated for appli- 
cation to a full-scale wing of different elasticity distribution. 


(2) It evaluates the variation of effective section lift-curve slope 
and of aerodynamic twist across the span as influenced by the presence of 
the fuselage, the nacelles, and other bodies on or near the wing (these 
obviously include spoilers, ailerons, and flaps). 


(5) It determines the manner in which section lift-curve slope varies 
with Mach number. 


Although the variation of lift-curve slope with Mach number can be 
obtained by the foregoing method by using wind-tunnel data at various 
Mach numbers, it is sometimes desired to determine the change (or changes) 
in lift distribution for section configurations for which suitable wind- 
tunnel data are not available, for example, variation of section 1ift 
with flap deflection Cis 


If the incompressible value of c7z can be obtained or estimated 


by some means (e.g., from tests for the flap deflected at some other sta- 
tion) and if the effective sweep angle at the new flap station is known, 
then compressible values of Clas can be obtained by substituting the 


Cla values for the section lift-curve slopes in equation (A359) 


m 
n= 


° " - M*cos?/w 


where My and m are the compressible and incompressible section lift- 
curve slopes, respectively. 
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It is therefore advantageous to evaluate the spanwise variation of 
Lage The effective sweep angle for each section can be obtained from 


equation (A39) by using the previously measured variation of my with 


Mach number. For example, equation (A39) can be rewritten as the linear 
equation (linear in squared terms 


Cs) + (io? ye )- 


This equation can be solved for m and Ay by a least-squares procedure 


if a series of equations are formed by substituting the values of Mmo 
and M obtained for each of several Mach numbers. 


Alternatively, a solution for m and Ay can be obtained by a 
graphical procedure (ref. 19). Plot the variation of mg with M in 
a fashion such that the abscissa (M-scale) is proportional to the fac- 


tor 1 / - MÉcos?Ay. Ihe plotted points will fall on a straight line 


which passes through the pole ms = O, 1A l- MPcos“Ay = o), provided 
Am is correctly chosen and the law (eq. (A39)) applies. Note that 


i/ "1 - MPcostAM = 1.0 when M= O. A form is presented in figure 9 by 


which the my values can be conveniently plotted on such an abscissa 
scale for each of a number of values of Aye For example, if a value 

of AM is tentatively chosen as Ay = 30°, draw a horizontal line inter- 
secting the right-hand ordinate at 30°. The intersections of this hori- 
zontal line. with the Mach number lines are the abscissa locations for 
the indicated discrete values of M. Values of mo are then appro- 
priately plotted vertically above or below these intersections. The 
intersection of this horizontal line with the pole-distance curve gives 
the pole location in scale units to the left of the abscissa point M - O. 
The desired value of Ay is then the one which gives the most nearly 
linear variation of the plotted. data with consideration given to the 
pole point. The desired m is the value of the intercept at M equal 
to zero. 


lliustrative example.- This section presents an analysis of wind- 
tunnel data obtained on a flexible model to eee the compressible 
section lift-curve slopes {m}. 


The model was a wing-fuselage configuration with nacelles mounted 
below and forward of the wing on ade struts attached to the wing 
semispan stations 7 = 0.37 and 0.65. The wing had the following addi- 
tional pertinent physical characteristics: 


ABDSCU FELIG s, ds e a A A AA a OOO 
Taper ratio + e s . * o . . e e > LÀ . *. . + LÀ LÀ . e a € o * . O. ho 
Sweep at quarter-chord dins. Ted aS “ee c^ e RR Koc a oe UO» www CD 


and the locus of aerodynamic centers was assumed to coincide with the 
quarter-chord line. 


——— -— — Rim — n ma pet ap e - - 


= —— —— —— —  )— a == —— A A O a es ds RS de a 
" : - "^ a eaves E" = re ge -3 * bum - 
EE ‘ f - = - a PU A ` 
$ . WU. a cma t aieo 24 Bc à PE sd ue tie v e 
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The model wing was instrumented with strain gages cemented to the 
surface of the wing steel spar along one semispan and with pressure 
orifices located along streamwise chord sections on the opposite semi- 
span at stations 7 = 0.155, 0.35, 0.56, 0.75, and 0.92. 


Tests were made at root-section angles of attack of 0° and 6° for 
Mach numbers 0.30, 0.50, 0.70, 0.75, and 0.80. 


The spanwise variation in section angle of attack for each test con- 
dition was obtained; from the root-section angles of attack in conjunction 
with strain-gage readings. The spanwise variation in section normal-force 
coefficient for each test condition was obtained by a spanwise fairing of 
the various local integrated chordwise pressure distributions. 


Elements of the [Sy | matrix were computed for wing semispan stations 


n = 0.10, 0.30, 0.50, 0.70, 0.85, 0.925, and 0.975. The resulting matrix 
S 

0.14869 -0.04630 -0.01249 -0.00482 -0.00149 -0.00082 -0.00063 
-0.04084 0.14811 -0.05547 -0.00793 -0.00188 -0.0009h -0.00070 
-0.00287 -0.01260 0.07827 -0.02713 -0.00285 -0.00121 -0. 00084 
-0.00042 -0.00068 -0.00426 0.04582 -0.01109 -0.00229 -0.00130 
-0.00013 -0.00017 -0.00053 -0.00522 0.04153 -0.01143 -0.00287 
-0.00008 -0.00009 -0.00023 -0.00093 -0.00617 0.03969 -0.01307 


-0.00006 -0.00006 -0.00015 -0.00050 -0.00124 -0.0079} 0.02917 


where the columns read down from the tip to the root and the rows read 
across from the tip to the root. For these same stations Ac,c/k and 
increments of section geometric angle of attack Aa due to the 6° change 
in root angle of attack were computed from data at each test Mach number 
and tabulated in table I. 


Substitution of {ar} = f^cnac for (2) and faat for {ee} in 


equation (12) gives 
O 
[S E25] - E (e) 


" 
Bata 
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O 
from which the values [m were computed and tabulated in table I. 
Values of m, for 7 = 0.975 are not shown since m. is very sensitive 


to the manner in which the spanwise variation of Ac, is faired in the 


region of the wing tip. These values for each wing station are plotted 
against Mach number in figure 10(a). 


Since the [94] matrix is based upon the "wing alone" configuration, 


the body interference effects which vary with body angle of attack will 
be apparent in the values of lift-curve slope particularly for those 
Stations near the wing root. 


lt was assumed that the Prandtl-Glauert relationship was applicable 
for these data. Accordingly the fmo} values were plotted (by using a 
form similar to that shown in fig. 9) for various selected values of Ay. 


The plots which yielded the most nearly linear variation of the data 
are reproduced in figure (10(b) for each span station. In selecting the 
plot which gave the most linear variation of mo, consideration was given 


to the fact that at very low Mach numbers (approx. 0.30) the value of mo 
could be affected by Reynolds number; whereas those obtained for high 
Mach numbers would reflect the effects of shocks, and so forth. These 
points were therefore given less weight in determining the best fit. 


Plots showing the spanwise variation of effective section sweep 
angle Ay and lift-curve slope m are given in figure 11. 
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APPENDIX H 
CALCULATION OF THE DIVERGENCE DYNAMIC PRESSURE 


Although the sweptback wing is usually considered to be divergence- 
free, it is conceivable that, with a large external diverging torque 
such as may be contributed by a tip tank, the wing could diverge. This 
appendix is concerned with the calculation of the dynamic pressure at 
which divergence will occur. 


An expression for the static lift on a flexible swept wing with an 
external tank has been derived in appendix D and, for a divergence inves- 
tigation, may be written as follows: 


E i se - [S2] 7 = {fe} m 


This is a matrix equation in which the elements in the column represent 
tbe lifts on the various spanwise segments of the wing. The total lift 
can be found from equation (HL) by multiplying the lift per unit span 


by the row matrix |2n] where 2h defines the segment width in the 
spanwise direction. Then 


L = 22m] (1) = alenl [E - Ja - [52| l {ox} 


or 


-e GE -e eq c 


Equation (H2) represents the lift on a wing in static equilibrium. 
Under the condition of wing divergence the equilibrium wing lift would 
be infinite. 


ee ee M 


|4H 


- —— ee M oe ee om occ Ttc n 
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The divergence Speed is then represented by the lowest value of q 
for which the lift in equation (H2) becomes infinite. For the right- 
hand side to be infinite, the determinant of the matrix whose inverse 
is given in equation (H2) must be zero, or 


EH i ARE - [82 ]| = o. 


which can be written in the form 


(Jel -aR e em 


The procedure is to solve for the lowest value of q which satisfies 
equation (H3). From the Cayley-Hamilton theorem, the dominant or highest 
modulus root A in the equation 


Nz] - [o] = [o] (8) 


may be found by iterating the matrix [D] . In this case equation (H3) 
may be put in the form of equation (H4) by multiplying through by 


Jl 


The result is 


E - a aaa 


— —M — — — —— 


—————————————— ~ TO————— n Á— e n 
+ 4 ^ ” . a * on « ^ 
4 - ` » = ry X- `a * - s - 
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Now A= D therefore, iteration of the matrix product 


Jel) jede] A es 


will yield the lowest value of q, which in this case is the desired 
divergence dynamic pressure. 


The theory of matrix iteration may be found in reference 20; however, 
a simple outline of the steps required to iterate the matrix product (H6) 
to obtain the divergence dynamic pressure is given here. First a trial 
column is chosen and this column is premultiplied by matrix (H6) to 
obtain a result column. The elements of this result column are divided 
by the last element of the result column and then become the elements of 
& Second trial column (the last element will be unity, having been 
divided by itself). The second trial column is then premultiplied by 
matrix (H6) to yield & second result column. This procedure is repeated 
until the same value is obtained for the last element in two successive 
result matrices. The reciprocal of this value is the desired divergence 


dynamic pressure q. 


-— — — AAA AR — - 
.- 
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TABLE I.- MEASURED DATA AND RESULTING SECTION LIFT-CURVE SLOPES FOR 


WIND-TUNNEL MODEL OF APPENDIX G 
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0.10125 


»10140 
. 10148 
. 10217 
.10317 
. 10h01 
¿10455 
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As 
_--Equilibrium position of 
section zero-lift line 
af 
ag 
ar Root section zero- 
q ift line 
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Figure 1.- Angle-of-attack definitions and sign conventions. Positive 
angles shown. 
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Figure 2.- Typical vortex locations and location of pertinent points on 


an arbitrary wing plan form. 
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Figure 3.- Structural skeleton of outboard sections of left wing shown 
in pl&n view. 
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Figure 4.- Plan view of left wing tip section. 
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Figure 5.- Polar representation of left wing used in development of the 
[So | matrix e 
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(Note that a straight-line image bound 


vortex is used instead of its true shape, 


(a) Front view of typical vortex-image system (high-midwing configuration). 


Figure 6.- Diagram of vortex-image system for & wing-fuselage combination. 
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(b) Plan view of typical vortex-image system (midwing configuration). 


Figure 6.- Concluded. 
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Figure 7.- Plot of overvelocity contours. 
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Figure 7.- Concluded. 
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Figure 8.- Diagram of a representative wing—tail-boom combination showing 
entry of tail load on to wing structure (d is distance between stations 
measured perpendicular to streamline). 
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Figure 9.- Plotting form based on equation m, = ————— for 


n - M*cos^Aw 
graphical determination of effective section sweep angle Ay for 
a known variation of lift-curve slope with Mach number. 
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Figure 10.- Variation at several semispan stations of section lift-curve 


slope with Mach number for linear scale and expanded M-scale. 
tunnel model of &ppendix G. 
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Figure 11.- Spanwise variation of effective sweep angle Ay and section 
lift-curve slope m for wind-tunnel model of &ppendix G. 
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